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I.  INTRODUCTION  TO  THE  RESEARCH  PROGRAM 


It  is  appropriate  to  begin  this  first  report  of  the 
series  with  an  overview  of  the  entire  research  program — its 
motivation,  objectives,  approach,  and  constraints.  This 
will  be  followed  by  a  summary  discussion  of  the  database, 
from  its  origin  to  its  production  on  digital  magnetic  tape 
and  in  hardcopy  display.  Finally,  detailed  measures  of  the 
database  occupy  the  bulk  of  this  report. 

STATEMENT  OF  THE  PROBLEM 

Recent  technological  developments  have  resulted  in  a  wide 
variety  of  imaging  systems  and  subsystems.  The  flexibility 
and  technologies  available  to  the  designer  include  various 
means  for  collection,  coding,  transmitting,  decoding,  analog 
and  digital  processing,  and  analog  and  digital  display.  The 
applications  of  such  systems  and  subsystems  are  myriad, 
tanging  from  static  and  dynamic  military  photointerpretive 
functions,  through  commercial  and  closed-circuit  television 
and  facsimile  systems,  to  diagnostic  radiological 
instrumentation  and  earth  resources  applications.  The 
scientific  world  is  quite  familiar  with  some  of  the 
techniques  which  can  be  used  to  "improve"  the  nature  of  any 
such  image,  and  the  non-scienti f ic  world  has  equally  seen 


such  as  the  Zapruder 


examples  of  processing  effectiveness, 
and  Hughes  films  of  the  Kennedy  assassination.  In  many 
cases,  it  is  clear  that  such  processing  and  display 

techniques  can  extract  information  in  the  original  image 
which  is  otherwise  well  below  the  threshold  capacity  of  the 
human  visual  system,  whereas  in  other  cases  it  is  clear  that 
processing  techniques  can  often  serve  either  to  hide 

existing,  and  important,  image  detail  or  to  "create"  image 
detail  which  is  perhaps  not  pr'  ^ent  in  the  original  image  or 

in  the  "real  world".  Heretofore,  most  of  these  areas  of 

image  system  and  subsystem  development  have  plainly  suffered 
from  their  inattention  to  human  observer  requirements.  This 
is  particularly  true  of  the  extensive  effort  in  digital 
image  processing,  especially  that  part  devoted  to  the 
improvement  ("enhancement",  "restoration")  of  images  for 
purposes  of  human  information  extraction.  In  nearly  all  of 
the  work  performed  in  laboratories  around  the  country  that 
are  pursuing  this  type  of  research,  the  necessary  evaluative 
efforts  to  determine  the  utility  of  processing  and  display 
techniques  have  not  been  conducted.  Rather,  reports  and 

publications  of  this  work  typically  take  the  form  of  "before 
and  after"  pairs  of  images,  where  the  reader  is  left  to 
estimate  the  utility  of  such  images  either  by  visual 
inspection  of  these  published  (second-  or  third-generation) 
photographs  or  by  the  subjective  opinions  offered  in  the 
text  by  the  author. 
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Because  the  intent  of  such  image  processing  techniques  is 
to  improve  the  information  extraction  capabilities  of  the 
human  observer,  it  is  clearly  appropriate  and  mandatory  that 
evaluative  techniques  include  objective  measurement  of  human 
information  extraction  from  such  images,  rather  than  merely 
subjective  estimates  of  the  overall  quality  or  utility  of 
the  image.  Unfortunately,  the  human  factors  experiments 
required  to  produce  quantitative  and  objective  assessment  of 
image  quality  have  rarely  been  conducted  in  image  processing 
laboratories  or  in  conjunction  with  image  processing 
programs . 

In  view  of  the  many  millions  of  dollars  being  devoted  to 
image  collection,  processing,  and  display  systems  for  the 
military  and  civilian  use  of  digitized  images,  it  is  quite 
clear  that  an  assessment  program  is  urgently  needed  to 
devise  procedures,  techniques,  and  metrics  of  digital  image 
quality.  Such  a  program  requires  the  establishment  of  a 
standardized  set  of  procedures  for  obtaining  human  observer 
information  extraction  performance;  relating  this 
performance,  in  a  quantitative  manner,  to  the  various 
collection,  processing,  and  display  techniques  and 
algorithms;  and  devising  a  quantitative  relationship  for  the 
multi-dimensional  scaling  of  the  various  collection, 
processing,  and  display  techniques  in  "performance  space". 

Only  through  such  an  integrated  program  of  research  can 
the  system  and  subsystem  designer  have  meaningful  data  for 


cost-benefit  analysis  of  future  system  development,  be  such 
systems  intended  either  for  military  or  for  non-military 
applications.  The  image  collection,  processing,  and  display 
technology  is  now  at  a  point  whereby  such  evaluative 
research  is  sorely  needed.  Fortunately,  microphotometric, 
microdensitometr ic ,  and  human  performance  measurement 
techniques  have  been  evolved  during  the  past  several  years 
to  relate  human  information  extraction  performance  to  the 
various  physical  characteristics  of  both  electro-optical  and 
photographic  image  displays.  The  present  research  program 
is  designed  to  extend  these  recently  developed  techniques 
into  the  arena  of  digital  images,  emphasizing  derivation  of 
metrics  of  image  quality  appropriate  to  digitized  images, 
and  providing  quantitative  cost-benefit  data  which  will 
permit  the  designer  and  system  developer  to  plan  his 
developmental  effort  as  well  as  to  specify  optimum  system 
components  for  particular  image  acquisition  and  display 
requirements. 


OVERVIEW  OF  THE  RESEARCH  PLAN 

The  research  plan  is  laid  out  schematically  in  Figure  1. 
Each  small,  solid-lined  box,  with  the  exception  of  the 
uppermost,  indicates  a  separate  task  to  be  conducted  during 
the  course  of  the  four-year  effort.  The  two  large,  broken- 
lined  boxes  delineate  the  specific  display  formats  that  will 
be  studied  and  compared  during  this  initial  program:  black 


•'  i 
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and  white  hard  copy  transparencies  and  electronic  displays, 
me  small,  broken-lined  box  at  the  bottom  illustrates 
important  extensions  of  this  research  to  be  pursued  in  the 
future,  namely  interactive  digital  displays  in  both  black 
and  white  and  in  full  color.  The  present  report  describes 
in  detail  the  execution  of  the  first  two  steps  of  this 
overall  plan,  specifically  the  creation  of  the  database  of 
known  content  and  with  controlled  image  degradation. 


PROGRAM  OVERVIEW 

Figure  1:  Schematic  diagram  of  proposed  research 

Research  Objectives 

The  research  objectives  of  this  program  are  as  follows: 

1.  Develop  standardized  procedures  and  techniques  to 
evaluate  hard-copy  (film)  and  soft-copy  (CRT)  digital 
image  quality. 
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2.  Compare  candidate  physical  metrics  of  image  quality. 

3.  Compare  hard-copy  with  soft-copy  displays  for  image 
interpretation. 

4.  Evaluate  candidate  processing,  enhancement,  and 
restoration  algorithms  for  improvement  of  image 
interpretation  on  soft-copy  displays. 

Specific  Research  Tasks 

In  keeping  with  the  general  goals  described  above,  the 
specific  research  tasks  are  as  follows: 

1.  Develop  an  imagery  database  and  image  interpretation 
scenarios  from  high  quality  aerial  photographs 
relevant  to  the  image  interpretation  task. 

2.  Select  and  purchase  display  and  interface  hardware  to 
present  the  image  database  on  soft-copy  (CRT) 
displays . 

3.  Develop  image  manipulation  software  for  soft-copy  and 
hard-copy  experiments. 

4.  Develop  and  standardize  observer  data  collection 
procedures  for  hard-copy  and  soft-copy  experiments. 

5.  Develop  and  standardize  procedures  for  obtaining 
physical  image  metrics  from  hard-copy  and  soft-copy 
displays . 

6.  Digitize  and  degrade  database  imagery  and  record 
images  on  hard  copy  and  magnetic  tapes  for  soft-copy 
display. 
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7.  Obtain  physical  image  metric  data  for  hard-copy  and 
soft-copy  displays. 

8.  Conduct  subjective  quality  scaling  and  information 
extraction  studies  on  hard-copy  images. 

9.  Conduct  subjective  quality  scaling  and  information 
extraction  studies  on  soft-copy  displays. 

10.  Evaluate  the  utility  of  image  quality  metrics  for 
both  hard-copy  and  soft-copy  imagery. 

11.  Conduct  subjective  quality  scaling  and  information 
extraction  studies  on  processed  soft-copy  imagery. 

12.  Analyze  utility  of  image  quality  metrics  for 

■* 

processed  soft-copy  imagery. 

13.  Compare  image  quality  metrics  for  hard-copy  and  soft- 
copy  (processed  and  nonpr ocessed )  images.  Relate 
these  results  to  concepts  and  models  of  human  visual 
performance  and  to  imaging  system  variables. 

Program  Ground  Rules 

The  overriding  objective  of  this  program  is  to  provide  a 
scientific  contribution  in  the  area  of  digital  image  quality 
and  image  metrics,  which  would  be  relevant  to  the  field  of 
military  photointerpretation  and,  at  the  same  time,  be 
available  to  the  scientific  community.  To  be  realistic  and 
consistent  with  tactical  interpretation  scenarios,  close 
coordination  was  maintained  with  Hq  USAF,  TAC,  and  RADC  at 
critical  stages  of  the  research,  in  particular  during 
imagery  selection. 
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The  degree  of  realism  was  also  enhanced  by  using 
photo  interpreters  with  experience  relevant  to  the  types  of 
imagery  employed.  Thus,  criteria  for  imagery  selection 
included  excellent  quality  and  content  typical  of  tactical 
scena  r  ios . 

Finally,  all  output  from  the  program,  including  results 
and  the  imagery  database,  are  unclassified  and  will  be 
reported  in  the  open  literature.  It  is  hoped  that  future 
investigations  of  digital  image  quality  will  use  the  same, 
available  database  to  permit  comparisons  with  and  extensions 
of  our  results. 
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II.  OVERVIEW  OF  DATABASE  DEVELOPMENT 

ORIGINAL  MATERIAL  AND  PHOTOINTERPRETATION  SCENARIOS 

The  original  negative  transparencies,  from  which  our  low- 
altitude  aerial  photographic  database  was  generated,  were 
acquired  from  the  film  libraries  of  the  Environmental 
Research  Institute  of  Michigan  (ERIM) .  With  the  support  and 
encouragement  of  sponsors  of  some  of  their  prior  research, 
Messrs.  Robert  Horvath  and  Stephen  Stewart  of  ERIM  aided  us 
in  locating  promising  candidate  materials.  These  were  taken 
to  Langley  AFB,  Virginia,  where  we  and  personnel  of  the 
460th  Reconnaissance  Technical  Squadron,  together  with  Dr. 
Harry  Snyder,  Dr.  Michael  Maddox,  and  Mr.  James  Turpin  of 
VPI&SU,  identified  suitable  7.6  x  7.6  cm  target  areas  and 
developed  essential  elements  of  information  (EEIs)  therein. 
Earlier  visits  to  Langley  AFB,  arranged  and  accompanied  by 
Mr.  Gilbert  Kuperman  of  AMRL,  Wr igh t-Pa tterson  AFB,  Ohio, 
had  resulted  in  the  full  and  enthusiastic  cooperation  of  the 
Tactical  Air  Command.  Their  support  and  encouragement 
continues  to  be  an  essential  contribution  to  the  viability 
of  our  studies. 

The  original  photographs  were  generated  by  a  K17  camera, 
with  a  74  degree  FOV,  W12  filter,  15.2  cm  foca.  length,  and 
1/250  shutter  speed,  at  f/11.  Four  are  in  Kodak  Plus-X 
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Aerial  Film  3401;  the  rest  are  in  Kodak  Tri-X  Aerographic 
Film  2403.  Of  the  10  scenes  used  in  our  hardcopy  studies,  4 
cover  U.S.  airfields,  4  more  cover  U.S.  naval  installations, 
and  2  illustrate  electronic  and/or  R  D  installations. 
Scales  range  from  1:2700  to  1:4400. 

Also  transported  to  Langley  AFB  on  several  occasions  were 
hardcopy  displays  of  digitally  processed  images,  degraded  by 
Gaussian  blur  and  noise.  These  were  prepared  by  Dr.  Robin 
Strickland  and  Mr.  Thomas  Bruegge,  employing  the  UA's  CYBEil 
175  computer,  the  DEC  11/70  -  IIS  70/E  pipeline  processor 
system  in  the  Digital  Image  Analysis  Laboratory  (DIAL),  and 
the  Eclipse  S/230-PDS  microdensitometer  in  the  Optical 
Sciences  Center  (OSC) .  The  purpose  of  these  early  blurred 
and  noisy  images  was  to  establish  a  meaningful  range  of 
these  parameters  for  use  in  subsequent  psychophysical 
studies.  It  was  our  intention  to  include  material  both 
better  and  worse  than  that  typically  encountered  by  tactical 
photointerpreters,  with  an  adequate  sampling  of  the 
parameter  space  between  these  extremes.  Several  iterations 
were  required  to  obtain  the  desired  range  and  sampling 
density . 

We  began  with  the  idea  of  using  image  contrast  as  a 
third,  independent,  dimension.  Discussions  with  senior 
interpreters  indicated,  however,  that  this  was  unnecessary, 
if  not  undesirable.  Something  more  akin  to  a  dimension 
measuring  signal-to-noise  ratio  (SNR)  was  indicated. 
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Because  all  scenes  contained  targets  of  varying  size  and 
contrast,  we  adopted  the  maximum  (peak)  intensity  as  our 
definition  of  "signal”,  adding  noise  proportional  to  peak 
signal.  All  scenes  were  scanned  to  yield  peak  digital 
values  of  approximately  2000.  Nominal  SNRs  of  200,  100,  50, 
25,  and  12.5  were  selected  for  these  studies. 


INITIAL  SELECTION  OF  EQUIPMENT  AND  MATERIALS  FOR 

DIGITIZATION 

Our  plans  called  for  the  OSC  Eclipse-S230/PDS- 
microdensitometer  to  serve  as  our  primary  facility  for  the 
generation  of  the  large-scale  (4096  x  4096)  digital  images 
on  magnetic  tape  (digitization  at  20  ym  steps  with  Gaussian 
scanning  apertures  of  20,  40,  80,  160,  and  320  urn  full-width 
half  maximum  (FWHM) ) ,  digital  noise  addition,  and  hardcopy 
display  (positive  transparencies  written  on  Kodak  Technical 
Pan  2415  by  the  microdensitometer  with  20  ym  aperture  and 
step  size).  The  DIAL  laboratory,  on  the  other  hand,  would 
serve  as  the  UA's  principal  facility  for  preparing  and 
displaying  softcopy  images,  in  512  x  512  format,  during 
future  studies  of  enhancement  and  restoration  techniques.  A 
similar  but  more  extensive  facility,  composed  of  a  DEC  11/55 
with  multiple  disk  drives  and  IIS  70/E  processor,  was 
developed  at  the  Human  Factors  Laboratory  at  VPI&SU  for 
softcopy  display  of  both  the  large  and  small  digital  image 
formats.  It  is  currently  in  use  for  soft-copy 

psychophysical  studies  paralleling  the  hardcopy  experiments 
to  be  described  in  a  subsequent  report. 
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By  the  end  of  summer,  1979,  a  dozen  interpretive 

scenarios  had  been  worked  out  at  Langley  AFB.  Ten  were 

selected  for  the  psychophysical  studies;  the  others  were  set 

aside  for  purposes  of  illustration  and  demonstration,  as  in 

this  report  (see  Appendix  B)  .  The  original  22.9  x  22.9  cm 

negatives  were  then  transported  to  Precision  Photo 

Laboratories,  Dayton,  Ohio,  where  two  high-quality  positive 

transparencies  of  each  were  prepared  by  contact  printing. 

During  the  printing,  each  image  was  gray  scale  stretched  to 

produce  approximately  the  same  maximum  dynamic  range  across 

all  scenes.  These  transparencies  became  our  "ground  truth" 

data.  Because  the  blurring  was  to  be  accomplished  by  analog 

means  (scanning  with  large  Gaussian  apertures  at  the  time  of 

digitization),  only  positive  transparencies  could  be  used 

without  incurring  prohibitive  computation  costs.  The 

positives  were  ready  at  the  end  of  September  1979; 

microdensi tometr ic  scanning  of  the  10  scenes  at  each  of  5 

o 

levels  of  blur  could  begin.  Each  of  the  50  (4096)  scans 
would  consume  3  hours.  An  additional  3  hours  would  be 
required  to  generate  each  of  the  250  (5  noise  levels)  first- 
generation  positive  transparencies  on  Kodak  Technical  Pan 
2415,  to  be  hand-processed  in  HC  110-D.  Five  second- 
generation  positives  of  each  of  the  250  transparencies  would 
then  be  prepared  on  Kodak  Duplicating  film,  SO-015,  by  hand 
processing  in  Dektol.  Hardcopy  psychophysical  studies  were 
scheduled  to  begin  in  mid-July,  1980,  and  had  to  be  finished 
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in  September.  We  needed  approximately  1000  hours  of 
faultless  performance  from  the  Eclipse/PDS  microdensitometer 
between  October  1979  and  July  1980.  It  was  not  to  be.  The 
digital  interface  (of  local  design  and  construction)  between 
the  microdensitometer  and  the  computer  proved  unreliable. 
We  had  completed  only  28  of  the  50  required  scans  and  none 
of  the  250  hardcopy  displays  when  we  had  to  abandon  it  in 
May  to  embark  on  a  "crash”  program  of  hardcopy  preparation. 

FINAL  PREPARATION  OF  DIGITAL  DATABASE  AND  HARDCOPY  DISPLAY 
Mr.  Ray  Schmidt  of  the  Image  Processing  Institute 
(University  of  Southern  California)  generously  offered  to 
take  on  the  task  of  generating  the  250  originals  from 
digital  tapes  supplied  by  us.  His  Dicomed  D-47,  DEC  11/40 
installation  was  already  prepared  to  routinely  write  4096  x 
4096  pixel  images  to  Kodak  2415  10.16  x  12.7  cm  sheet  film 
in  86  mm  format,  and  he  had  extensive  experience  in  hand¬ 
processing  this  film  in  Kodak  D-19  developer.  Though  the 
resolution  of  the  Dicomed  is  such  as  to  provide  no  more  than 
3000  x  3000  independent  pixels,  our  time  constraints  were 
such  as  to  encourage  acceptance  of  this  limitation.  A 
second  major  problem  attendant  to  our  equipment  failure  was 
solved  by  Mr.  H.  D.  Fisher  of  the  UA's  Health  Sciences 
Center.  The  analog  blurring  operation  via  the 
microdensitometer  was  now  impossible.  Mr.  Fisher  offered  to 
do  the  blurring  digitally  by  overlapping  9  x  9  (512) 2  fast 
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Fourier  transforms,  each  done  "in  place"  in  the  large  memory 
of  his  VAX  11/780.  Each  of  the  10  (4096) 2  original 
microdensitometer  scans  of  the  Precision  Photo  positives 
would  be  read  to  disk,  Fourier  transformed,  and  multiplied 
by  four  appropriate  Gaussian  filter  functions,  thereby 
creating  the  spectral  versions  of  the  blurred  images.  These 
would  be  inverse  transformed  to  yield,  after  discarding 
overlap,  the  required  images.  With  the  very  efficient 
software  written  by  Messrs.  Fisher  and  Kevin  McNeill,  each 
scene  could  be  so  processed  in  four  hours  by  the  VAX.  The 
blurred  image  files  were  then  written  to  tape  in  internal 
format,  one  16-bit  word  per  pixel,  for  subsequent  input  to 
programs  that  affected  noise  additions  (five  levels),  look¬ 
up  table  transformations,  and  truncation  to  8-bit  pixels. 
Thus  were  the  digital  tapes  prepared  for  display  at  the 
Image  Processing  Institute. 

The  Dicomed-wr itten  originals  were  done  in  "log 
complement"  mode  to  yield  the  positives  required  by  our 
duplicating  process.  The  selection  of  "log"  mode  for 
Dicomed  film  writing  was  consistent  with  the  experience  of 
the  Image  Processing  Institute.  The  look-up  tables  were 
designed  in  such  a  way  as  to  reproduce  the  contrast  of  the 
original  Precision  Photo  positives  in  the  case  of 
uncorrupted  images  (10  in  number),  allowing  the  remaining  24 
versions  of  each  scene  to  take  on  the  contrast  reductions 
attendant  to  both  blurring  and  noise  addition.  The  final 


copies  fell  somewhat  short  of  this  design  goal  because  of 
inconsistencies  in  the  exposure  and  development  of  both  the 
250  originals  and  the  1250  copies.  As  will  be  discussed  in 
a  subsequent  report,  these  inconsistencies  do  not  appear  to 
have  caused  any  first-order  effects  in  the  results  of  the 
psychophysical  studies. 

SUMMARY  MEASURES  OF  DIGITAL  DATABASE  AND  HARDCOPY  DISPLAY 
In  the  following  two  tables  we  list  the  nominal  and 
actual  values  of  blur  and  SNR  (peak  signal  divided  by  rms 
noise)  exhibited  by  the  digital  data  recorded  on  magnetic 
tape,  on  the  one  hand,  and  its  display  on  SO-015,  on  the 
other.  The  differences  result  from  undesired  blurring  and 
noise  addition  by  the  Dicomed  film  recorder. 


TABLE  1 

Nominal  and  Actual  Blur  Values 


Level 

Blur  Dimension 

Nominal 

(FWHM  in  pm) 

Tape 

SO  015 

0 

20 

22 

40 

1 

40 

43 

52 

2 

JO 

81 

84 

3 

160 

161 

162 

4 

320 

320 

322 
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TABLE  2 

Nominal  and  Actual  Signal-to-No ise  Levels 


Level 

Noise  Dimension 

Nominal 

Tape 

SO  015 

0 

200 

200 

75 

1 

100 

100 

60 

2 

50 

50 

42 

3 

25 

25 

24 

4 

12.5 

12.5 

12 

These  values  of  additive  noise  do  not  reflect  the 
contribution  of  the  {signal  dependent)  grain  noise  exhibited 
by  the  positive  Precision  Photo  transparencies,  which 
represent  our  "ground-truth"  data.  This  contribution  is 
discussed  later  in  the  text. 

AVAILABILITY  OF  DATABASE 

The  digital  database,  in  the  form  of  binary  magnetic 
tapes  with  4096  bytes/r ecord ,  1  byte/pixel,  and  4096 
records/image  will  be  available  to  other  investigators  when 
this  program  terminates.  We  will  not  display  it  in  this  or 
subsequent  reports,  so  that  future  subjects  of 
psychophysical  experiments  will  have  had  no  prior  exposure 
to  it.  For  the  same  reason,  we  do  not  discuss  its  content, 
except  in  most  general  terms. 
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The  reader  may  gain  an  impression  of  the  range  of  blur 
and  noise  covered  in  the  database  by  referring  to  Appendix 
B,  where  we  show  degraded  versions  of  an  aerial  photograph 
of  the  US  Army  Tank  and  Automotive  Command  (TACOM) ,  Warren, 
Michigan. 
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III.  SELECTION  OF  QUALITY  DIMENSIONS  OF  DATABASE 

As  noted  in  the  preceding  overview,  a  relevant  range  of 
quality  parameters  was  established  in  several  iterations, 
each  involving  the  preparation  of  a  small  set  of  digitally 
degraded  images,  followed  by  presentation  to 
photointerpreters  (Pis)  in  order  to  determine  the  maximum, 
minimum,  and  intermediate  values  of  resolution,  noise,  and 
contrast.  It  was  our  intention  at  the  outset  of  the  program 
to  consider  six  levels  of  blur,  three  levels  of  noise,  and 
three  levels  of  contrast.  Each  of  these  quality  parameters 
independently  affects  ’•signal"  as  perceived  in  the  image. 
Feedback  from  Pis  suggested  that  two  dimensions  of  quality 
reflecting  resolution  and  signal-to-noise  ratio  would  be 
adequate  and  realistic.  Consequently,  image  contrast  was 
rejected  as  an  independent  dimension,  since  Pis  typically 


work  with 

transparencies  or 

electronic 

d i spl ays 

which  are 

optimally 

contrast-stretched 

for  viewing 

.  A  final  decision 

was  made 

to  adopt  five  levels  of  blur 

and  five 

levels  of 

noise . 
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Pis  indicated  that  an  image  quality  range  of  4  to  8  (or 
9)  on  the  0-9  NATO  rating  scale  is  typical  of  their  working 
environment.  In  terms  of  ground  resolved  distance  (GRD) , 
this  is  a  spread  of  approximately  1.2-1. 5  m  to  5-15  cm, 
respectively  (see  Appendix  A).  Estimating  our  ground  truth 
GRD  in  the  Precision  Photo  positive  transparencies  to  be  in 
the  region  of  5-8  cm,  we  adopted  five  degrees  of  blur, 
giving  GRD  values  of  approximately  7  cm,  15  cm,  30  cm,  0.6 
m,  and  1.2  m.  The  factor  of  two  between  blur  levels  agrees 
with  the  NATO  and  other  rating  scales,  which  correlate  with 
logarithmic  GRD  increments. 

NOISE 

Consultation  with  Pis  established  a  maximum  rms  level  of 
Gaussian  additive  noise  corresponding  to  a  lower  limit  to 
our  signal-to-noise  ratio  dimension  of  12.5,  where  signal- 
to-noise  is  defined  as: 

SNR  =  (peak  signal)/(rms  noise).  (1) 

The  selection  of  intermediate  levels  in  the  noise  dimension 
was  based  on  the  desire  to  create  subjectively  equal 
increments  in  noise.  It  is  consistent  with  Weber's  law  of 
psychophysical  theory  that  subjectively  equal  steps  in 
random  grain  noise  should  occur  at  logarithmic  intervals. 
Zwick  and  Brothers  (1977)  report  that  just-noticeable 
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differences  in  rms  granularity  occur  for  intervals  of  6,  16, 
or  30%,  depending  on  whether  the  scene  is  uniform, 
"average,"  or  "complex",  respectively.  Tests  involving  the 
addition  of  Gaussian,  uncorrelated  noise  at  different 
increments  in  rms  level,  a  ,  indicated  that  a  factor  of  two 
is  suitable  for  coverage  of  the  noise  dimension; 
accordingly,  five  values  of  SNR  were  selected:  ground  truth 
(=  200),  100,  50,  25,  and  12.5. 


CONTRAST 

The  positive  transparencies  produced  from  the  original 
ERIM  negatives  by  Precision  Photo  Laboratories  are  of  high 
quality  in  terms  of  their  tone  rendition.  They  had  been 
prepared  in  such  a  way  as  to  satisfy  the  familiar  PI 
criterion:  "Stretch  the  contrast  as  far  as  you  can  without 
losing  anything."  We  experimented  with  several  digital 
algorithms,  particularly  histogram  equalization,  as 
potential  candidates  for  a  universal,  scene-independent 
preprocessor.  None  was  found  that  did  not  noticeably 
suppress  important  details  in  one  or  more  scenes.  We 
accordingly  decided  to  use  the  inherent  contrast  of  the 
Precision  Photo  transparencies  as  the  design  goal  of  our 
displays.  The  undegraded  digital  images  would  thus  be 
displayed  to  match  the  "ground  truth"  data  in  the  contrast 
dimension.  The  degraded  images  would  then  manifest  the 
contrast  reductions  inherent  to  the  operations  of  blurring 
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and  adding  noise.  Scene-dependent  contrast  stretching,  on 
the  other  hand,  would  be  regarded  as  a  form  of  digital 
processing,  to  be  studied  in  another  phase  of  the  program. 
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IV.  GENERATING  THE  DIGITAL  DATABASE  FROM 
ORIGINAL  POSITIVE  TRANSPARENCIES 

FORMAT  OF  THE  DIGITAL  DATABASE 

The  area  scanned  from  the  Precision  Photo  transparencies 
for  each  scene  was  82  x  82  mm,  this  being  approximately  the 
size  necessary  to  be  representative  of  the  scene  content 
found  in  tactical  photo  interpretation.  The  sampling 
aperture  and  interval  for  digitization  was  20  pm,  yielding  a 
"ground  truth"  resolution  on  the  order  of  7.6  cm.  Each 
original  scene  thus  provided  a  matrix  of  4096  x  4096  pixels. 

Each  pixel  was  ultimately  scaled  down  to  8  bits  or  256 
intensity  levels  after  the  digitization  and  computer  data 
processing  described  later.  This  number  of  quantization 
levels,  besides  being  adequate  for  psychophysical  studies, 
is  the  standard  in  most  digital  image  processing  hardware, 
for  example,  the  IIS  70/E  at  the  Digital  Image  Analysis 
Laboratory . 

A  stepwedge  was  combined  with  each  image  at  the  final 
stage  of  blurring  and  noise-addition  for  the  purpose  of 
photographic  calibration  in  hardcopy  production.  The 
stepwedge  replaces  the  last  128  lines  of  each  scene, 
corresponding  to  a  strip  of  2.5  mm  width  at  the  edge  of  the 
hardcopy  format;  this  strip  is  divided  into  16  blocks  having 
input  intensity  levels  0,  17,  34,  51,  ...,  255. 
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The  block  diagram  of  Figure  2  outlines  the  steps  of  the 
database  production,  the  elements  of  which  are  discussed 


below. 
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TRANS¬ 
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DIGITIZE 

4+096x^096 
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INTERVAL 
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10  IMAGES 
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250  IMAGES 
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RECORD 

1(096  RECORDS 


FFT  FILTERING  NOISE  ADDITION  AND  SCALING 


(n  pm  =  GAUSSIAN 
FWHM) 


TO  8  BITS 


VAX  11/780 

Figure  2:  Digital  database  production 


IMAGE  IDENTIFICATION  CODE,  ORDER-OF-BATTLE 

A  four-digit  code  identifies  each  degraded  image  as 
follows : 


1st  Digit:  Photo 


0,1, 4, 5, 6, 8, 9 


2nd  Digit:  Area  within  photo  1  or  2 


3rd  Digit:  Blur  level 


0-4,  increasing  in  blur 


4th  Digit:  Noise  level 


0-4,  increasing  in  noise 
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Twc  non-overlapping  areas  were  selected  from  each  of 
three  Precision  Photo  transparencies  identified  as  photos  1, 
4,  and  5.  The  order  of  battle  for  each  photo  is: 

Air:  0 ,  1 ,  6 

Sea:  4,  5 

Land:  8 ,  9 

Throughout  this  report  we  often  refer  to  scene  41,  say, 
as  a  convenient  abbreviation  of  photo  4 ,  area  1 . 

SCAN  CHARACTERISTICS  OF  THE  PERKIN-EI MER  PDS 
*!ICR3DENSIT0NETER. 

The  PDS  is  capable  of  generating  12  bits/pixel,  i.e., 
409S  gray  levels.  An  11-bit  dynamic  range,  yielding  values 
from  0  to  2000  was  convenient  for  our  d ig i t i za t ion .  The  PDS 
was  operated  in  the  transmission  mode  (number  proportional 
to  transmission),  with  the  highest  transmission  region  or 
highlight  within  the  scene  calibrated  to  N  =  2000.  Figure  3 
shows  characteristic  curves  obtained  by  scanning  a  Kodak 
density  wedge  with  different  settings  for  the  minimum 
density.  The  gain  of  the  characteristics  is  essentially 
unchanged  as  the  highlight  calibration  varies;  it  is 
described  by: 

D  ~  ^Dmin  “  Dmax^lcg  Nmax  ^  lc^l0N  +  Dmax  * 
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The  above  expression  is  valid  over  a  range  of  log  N  of  2.3 


(as  specified  by  Pe r kin-Elmer ) .  Appendix  C  describes  in 
more  detail  our  use  of  the  microdensitometer  and  explains 
the  nonlinearities  in  Figure  3. 


ORIGINAL  FILM  DENSITY 

Figure  3:  PDS  microdensitometer  scan  characteristic 

Ground  truth  images  of  4096  x  4096  pixels  were  scanned  at 
20  um  pixel  intervals  using  a  square  area  of  dimension  8.19 
cm  on  the  Precision  Photo  original  transparencies.  The 
resolution  of  this  scanning  aperture  is  shown  in  the  MTF 
calibration  curve  of  Figure  4,  obtained  by  scanning  bar 
patterns  of  known  width.  Note  that  the  minimum  of  this 
approximate  sin(x)/x  distribution  occurs  at  a  spatial 
frequency  of  45  cycles  mm"^,  This  measurement  reflects  the 
combined  effects  of  the  size  of  the  physical  aperture  and 


/ 


the  image  degradation  associated  with  the  optics  of  the 
microdensitometer . 


LINES/mm 

Figure  4;  PDS  MTF  (20  pm  aperture) 

GRAIN  NOISE  IN  PRECISION  PHOTO  TRANSPARENCIES 

The  10  original  scenes  were  scanned  from  positive 
transparencies  made  by  contact-printing  7  ERIM  negative 
transparencies,  the  film  type  being  either  PLUS-X  ?r  TRI-X, 
as  shown  in  Table  3. 

Film  grain  noise  in  the  positives  varies  according  to: 

1.  The  granularity  of  the  negative  emulsion  (TRI-X  is 
rated  1.65  times  as  grainy  as  PLUS-X),  and 

2.  The  degree  to  which  the  negatives  (of  varying 
contrasts  and  mean  density)  were  contrast-stretched 
during  printing  to  achieve  positives  of  equal 
contrast. 
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TABLE  3 

Film  Type  for  Each  Image 

ERIM  Negative  Film  Photo  #,  Area  # 

01 

11 

PLUS-X 

12 

61 

41 

42 

51 

TRI-X 

52 
81 
91 

The  classical  model  of  film-grain  noise  predicts  a 
square-root  law  between  the  rms  density,  oD,  and  the  mean 
density,  D 
thus : 

°D  =  k(D)1/2.  (3) 

The  constant  k  is  the  value  of  granularity  at  unit  density. 
The  resulting  noise  in  the  corresponding  digital  value,  N, 
is  derived  using: 


27 


log10  N  -  1Dmax  -  TD  <4> 

°N  =  CdN/dDJo  D.  (5) 

Thus 

loNl  =[YNk/0.43433L(Dmax  -  ( 1  ogj  QU)  /  y)  )]1/2 ,  (6) 

or, 

N/loNl  =  [0. 4343/Yk]  l(Dmax  -  (log  1()N/y)]~1/2,  (7) 

where : 

0  <  N  <  255 

Y  =  gamma  of  PDS  characteristic. 

Figure  5  gives  curves  of  " signal- to-noise"  (N/oN)  versus 
N  for  each  photograph.  These  were  obtained  by  varying  k  in 
order  to  test  the  fit  of  the  expression  for  N/|oN|  to  three 
or  four  measurements  of  oM  taken  at  int'ervals  in  N  between  0 
and  255  in  local  image  regions  of  apparently  uniform 
intensity  (i.e.,  no  visible  structure  apart  from  grain 
noise . ) 

In  aodition  to  the  intentional  noise  of  four  levels  (o„  = 
20,  40,  80,  160),  we  therefore  have  a  source  of  signal- 

dependent  noise  to  consider  in  the  digital  database.  The 
fact  that  this  noise  is  introduced  before  blurring  means 
that  it  is  halved  with  each  increasing  level  of  blur.  A 
third,  additive  source  of  noise  in  the  hardcopy  database  is 
discussed  in  Appendix  E. 
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Figure  5:  Grain  noise  in  Precision  Photo  transparencies,  o, 


=  rms  noise 


COMPUTER  PRODUCTION  OF  BLURRED,  NOISY  IMAGES  FROM 
GROUNDTRUTH  BLUR 


Blur 

It  was  our  original  intention  to  accomplish  blurring  by 
scanning  the  Precision  Photo  transparencies  using  Gaussian 
apertures  manufactured  by  exposing  the  necessary 
transmission  distributions  on  film.  The  four  blur-widths 
corresponding  to  our  range  of  GRD  would  be:  R^  =  40,  80, 
1*0,  320  ^m,  where  R^  is  the  full-width  at  half  maximum 
{ FWHM) .  The  usual  measure  of  Gaussian  width  ob  (standard 
deviation)  is  related  to  Rh  by  *  0.425  R^ .  Because  the 
images  are  sampled  on  the  PDS  at  20  urn  intervals,  it  is 
useful  to  specify  values  of  o.  as  0.849,  1.699,  3 .  3 9 8 r  and 
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6.795,  measured  in  units  of  inter-pixel  distance.  In 
reality,  analog  blurring  proved  impossible  due  to  equipment 
failure  causing  us  to  adopt  a  digital  method. 

The  technique  of  block-mode  filtering  by  Fast  Fourier 
Transforms  was  implemented  by  H.  D.  Fisher  of  the  University 
of  Arizona's  Health  Sciences  Center  on  the  VAX  11/780 
computer.  This  method,  shown  in  schematic  form  in  Figure  6, 
is  equivalent  to  convolving  the  original  digitized  data  with 
four  64  x  64  Gaussian  point  spread  functions  of  widths  ^ 
given  above.  In  block-mode  filtering,  each  of  81  (i.e.,  9  x 
9)  partially  overlapping  blocks  of  512  x  512  pixels  was 
Fourier  transformed  in  the  VAX  memory  and  multiplied  by  the 
appropriate  filter  (the  Fourier  transform  of  the  aussian 
point  spread  function  in  our  case)  .  After  inverse 
transformation,  the  32-pixel  wide  border  is  discarded  and 
output  blocks  of  448  x  448  pixels  are  set  out  contiguously 
as  in  Figure  6.  The  exception  to  this  occurs  at  the  edges 
of  the  4096  x  4096  matrix,  where  only  the  inside  borders  of 
the  512  x  512  blocks  are  discarded,  leaving  output  blocks  i  i 
480  x  480  pixels  in  each  corner  and  480  x  448  pixels  along 
the  edges.  The  expected  wraparound  effect  in  a  32-pixel 
wide  strip  at  the  outermost  border  is  not  evident  in  our 
images.  The  computer  processing  time  expended  in  generating 
the  4  blurred  versions  of  each  nc  the  10  409b  x  4096  scenes 
was  4  hours. 
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Figure  6:  Blurring  by  block-mode  filtering 


We  have  determined  by  experiment  that  the  resulting 
effective  GRD  is  related  to  ob  by 

GRD  =  2.5  u'b  (3) 

This  relation  is  valid  even  for  truncated  Gaussian  point 
spread  functions,  provided  that  the  extent  of  the  function 
is  at  least  3  oh>  For  the  case  of  airport  #1,  Scene  01,  for 


example,  we  therefore  have  GRDs  of  ground  truth,  18,  35,  71, 
and  142  cm  (based  on  a  photographic  scale  factor  of  4400:1 
for  this  scene)  . 


v  c  i  s  e . 


The  next  stage  in  image  g-  •  at  ion  u-.i 
generation  cf  a  single  digital  file  of  4091  x  40 91  pixels  cl 
random  Gaussian  noise,  which  was  then  to  be  weighted  -m  1 
added  to  each  blurred  scene.  During  initial  tests  with  12 
x  512  image  arrays,  we  created  pairs  of  uniformly 
distributed  random  numbers  (rv  and  rv,,)  using  the  IBM 
FORTRAN  routine  RANDIJ  for  subsequent  computation  of  Gaussian 
random  numbers  via  the  following  algorithm: 

g  =  a  (-2  log^  rv^)  cos  (2  n  rv.,)  .  (9) 

It  was  found,  however,  that  this  algorithm  periodically 
venerates  highly  correlated  structure  ever  an  area  of  4091  x 
4091  pixels.  (Forsythe,  Malcolm,  and  Mclcr,  1977,  report  a 
similar  phenomenon,  and  have  written  a  superior  algorithm 
'JR AND.  )  We  subsequently  generated  the  large  uncorrelated 
noise  file  more  economically  by  scanning  a  uniformly  exposed 
photographic  emulsion,  developed  to  a  density  cf  unity, 
using  a  small  (10  urn)  aperture  in  the  PDS  nicrcdensitcmeter  . 
After  subtracting  the  mean  value,  the  noise  statistics 
remain  approximately  Gaussian  (see  Figure  7),  where  the 
standard  deviation  is  approximately  40  intensity  units, 
’donee,  with  the  original  database  still  within  the  range  0 
tc  2000,  noise  weighting  factors  cf  0.5,  1,  2,  and  4, 
aopl ied  to  this  noise  file,  provided  the  desired  values  of 
additive  noise. 
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Figure  7:  Noise  histogram 

SUMMARY  MEASUREMENTS  OF  QUALITY  AND  STATISTICS  OF  DIGITAL 
DATABASE 

As  previously  mentioned,  the  digital  database  is  in  8-bit 
format  on  9-track  tape  with  4096  bytes  (pixels) /record  and 
4096  records/ image  file. 

Blur 

Table  4  compares  actual  and  nominal  values  for  the 
effective  FWHM  of  the  Gaussian  point  spread  functions  for 
each  level  of  blur.  The  small  deviations  from  nominal  are 
caused  by  the  combined  effects  of  the  PDS  scanning  aperture 
and  the  imaging  optics. 
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TABLE  4 


Nominal  and  Actual  FWHM  (urn) 


Level  Nominal  Blur  Digital  Database  on  Magnetic  Tape 


0  20  22 

1  40  43 

2  80  81 

3  160  160 

4  320  320 


Noise 

The  level  of  noise  introduced  during  scanning  by  the  PDS 
microdensitometer  is  negligible.  The  four  nominal  values 
for  rms  noise  were  obtained:  on  =  2.5,  5,  10,  20  (8-bit 
scale)  . 

Histogram-Related  Statistics 

We  used  a  512  x  512  subsampling  of  each  of  the  4096  x 
4096  images  to  derive  global  histogram  statistics.  Figure  8 
and  Tables  5  and  6  show  maxima,  minima,  and  means  for  blur 
levels  0,  2,  and  4.  They  also  list  the  values  of  a 
parameter  called  %-dynamic  range,  this  being  a  measure  of 
the  limiting  gray  levels  above  and  below  which  reside  x%  of 
the  total  number  of  pixels.  For  example,  90%-dynamic  range 
at  194  and  14  means  that  5%  of  all  pixels  have  gray  levels 
less  than  14  and  5%  have  gray  levels  greater  than  194. 
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Fiqure  8:  Illustration  of  %  D.R.  (dynamic*  range) 

Image  histograms  for  the  10  scenes  at  blur  level  0  are 
shown  in  Figure  9,  and  histograms  for  scenes  01,  41,  and  91 
at  blur  levels  2  and  4  are  given  in  Figure  10  (a-f) . 
Several  histograms  are  noticeably  bi-  or  tri-modal;  all  are 
biased  toward  the  zero  signal  or  high  density  end  of  the 
scale.  The  port  scenes  from  photo  4  (scenes  41  and  42) 
exhibit  somewhat  broader  gray-level  distributions, 
particularly  scene  41.  Figure  11  illustrates  the  potential 
utility  of  the  %-dynamic  range  values  as  measures  of  the 
image  "contrast"  or  "peak-to-peak"  signal.  For  example, 
scenes  81,  91,  and  12  possess  almost  identical  true  peak-to- 
peak  values,  yet  the  spread  in  their  90%  dynamic  range 
figures  is  large.  The  higher  proportion  of  highlight 


densities  for  scene  12  is  also  evident  in  the  shape  of  the 


TABLE  5 


First  Order  Statistics  for  10  Scenes,  Blur  0,  Noise  0. 


Image  # 

Max; 

Min 

Mean 

98% 

D.R. 

90% 

D.R. 

80% 

D.R. 

0100 

233 

0 

38 

158 

2 

103 

2 

89 

2 

1100 

227 

0 

23 

122 

2 

91 

3 

67 

4 

4100 

255 

0 

51 

140 

7 

109 

9 

96 

12 

5100 

207 

0 

17 

131 

2 

82 

3 

45 

3 

6100 

255 

0 

70 

255 

1 

235 

2 

176 

14 

8100 

244 

0 

11 

99 

2 

37 

3 

21 

3 

9100 

240 

0 

28 

173 

2 

92 

3 

72 

4 

1200 

239 

0 

30 

179 

2 

109 

3 

88 

4 

4200 

210 

0 

28 

93 

5 

72 

6 

61 

7 

5200 

203 

0 

11 

83 

2 

45 

2 

31 

3 
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TABLE  6 


Histogram  Statistics  for  Representative  Scenes,  Blurs  2 
(top)  and  4  (bottom)  ,  Noise  0. 


Image  t 

Max; 

Min 

Mean 

98% 

D.R. 

90% 

D.R. 

30% 

D.R. 

0100 

228 

0 

28 

156 

2 

101 

2 

88 

2 

4120 

255 

4 

51 

125 

7 

102 

9 

92 

12 

9120 

233 

1 

28 

170 

3 

90 

3 

71 

4 

0140 

244 

0 

38 

150 

2 

99 

2 

86 

O 

c. 

4140 

255 

4 

51 

119 

8 

100 

10 

91 

13 

9140 

230 

1 

28 

162 

3 

87 

3 

69 

5 

histogram.  These  histograms  illustrate  the  familiar 
difficulties  that  arise  when  one  attempts  to  find  a 
universal,  psychcphysically  meaningful  definition  of 
"signal"  in  an  optical  image. 

We  recall  from  Figure  2  that  there  is  no  renormalization 
after  the  blurring  process  ( FFT- f i 1 t er ing) ;  hence,  we  would 
expect  a  post-convolution  histogram  to  be  slightly  reduced 
in  dynamic  range,  and  its  shape  altered  somewhat, 
particularly  in  the  highlight  range  where  structure  tends  to 
be  fine  and  therefore  most  susceptible  to  blurring.  Figure 
12  illustrates  the  effect  of  blur  (levels  2  and  4)  on  the 
dynamic  ranges  for  scenes  01,  41,  and  91.  There  appears  to 
be  no  major  scene-dependent  effect  caused  by  blurring  of 
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Figure  10: 


Representative  digital  database  histograms 
{Blurs  2,  4:  Noise  0) 
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Figure  11:  Peak-to-peak  and  90%  dynamic  range.  Blur  0 

highlights,  since  the  spread  of  dynamic  range  with  blur  for 
100,  98,  and  90%  levels  is  approximately  constant  within 

scenes.  All  three  measures  of  dynamic  range  remain  roughly 
constant  throughout  the  blur  dimension.  We  do  not  reproduce 
here  the  histogram  for  blurs  2  and  4  for  every  scene  because 
the  relationship  among  their  dynamic  range  values  with  blur 
is  essentially  equivalent  to  those  shown. 

In  understanding  the  shape  of  a  histogram  after  the ^ 
addition  of  noise  there  are  two  factors  to  consider.  First, 
the  histogram  of  the  combined  image  (image  file  plus  noise 
file)  is  akin  to  the  convolution  of  the  histograms  of  the 
separate  components.  The  simplest  example  is  the  addition 
of  zero-mean  Gaussian  noise  to  an  image  of  uniform 
intensity,  k;  the  composite  histogram  is  obtained  by 
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Figure  12:  Effect  cf  blur  cn  dynamic  range,  zero  noise 


convolving  a  Gaussian  function  with  a  delta  function,  the 
result  being  a  Gaussian  histogram  shifted  to  r.  mean  level  k. 
Notice  also  that  the  inherent  smoothing  associated  with 
convolution  causes  considerable  reduction  of  the  absolute 
value  of  the  peak  in  the  uniform  histogram. 

The  second  factor  results  from  the  fact  that  noise 
addition  was  performed  by  combining  the  11-bit  (positive) 
signal  data  with  11-bit  (positive  and  negative)  noise  data. 
The  resulting  data  were  lmited  in  dynamic  range  by  setting 
negative  values  equal  to  0  and  values  exceeding  4095  (12 
bits)  equal  to  4095.  This  "scaling"  accentuates  the 
histogram  at  the  end  of  the  gray  level  range  towards  which 
the  original  histogram  was  biased.  Both  factors  are 


apparent  in  the  histogram  of  scene  91  at  blur  level  0  and 


noise  level  4,  shown  in  Figure  13.  When  scaled  to  8  bits,  a 
total  of  24%  of  the  image  pixels  are  at  level  0  due  to  the 
addition  of  noise  of  standard  deviation  20  (3-bit  scale)  to 
an  image  whose  mean  level  is  28. 


Figure  13:  Digital  database  histogram  (Noise  4) 
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V.  GENERATING  THE  HARDCOPY  DATABASE 

FORMAT  OF  HARDCOPY 

The  task  of  generating  hardcopy  was  transferred  from  the 
UA  PDS  microdensitometer  to  the  Dicomed  D-47  at  the  USC 
Image  Processing  Institute.  Our  specifications  called  for 
hardcopy  display  of  positive  transpa renc ies  written  on  Kodak 
Technical  Pan  2415  at  the  same  scale  as  the  original 
Precision  Photo  originals.  Playback  should  therefore  have 
been  at  pixel  spacing  of  20  urn,  equal  to  PDS  digitization. 
In  reality,  the  Dicomed  operates  at  21  wm  pixel  spacing, 
yielding  slightly  larger  hardcopy,  having  a  final 
reproduction  scale  of  36  mm:82  mm  or  1.05:1.  A  requirement 
for  five  copies  of  each  2415  transparency  for  PI  experiments 
necessitated  positive  contact  printing  on  Kodak  S0-015 
duplicating  film.  The  format  of  the  hardcopy  films  is  shown 
in  Figure  14. 


PLAYBACK  CHARACTERISTICS  OF  THE  DICOMED 

The  Dicomed  D-47  image  writer,  operating  in  its  "high 
resolution"  mode,  generates  a  matrix  of  pixels  spaced  at 
intervals  of  21  ;;m,  each  pixel  having  one  of  256  gray  levels 
corresponding  to  an  ultimate  film  density,  in  Kodak 
Technical  Pan  2415,  determined  by  various  electronic 
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Figure  14:  Hardcopy  format,  Kodak  2415 

calibrations  of  the  instrument  as  well  as  the  final  darkroom 
processing  in  D-19  developer.  Appendix  D  gives  the 

manufacturer's  specifications  for  this  instrument. 

Resolution,  MTF. 

Dicomed  resolution  is  specified  by  the  manufacturer  to  be 
3000  lines  over  the  plotting  matrix,  equivalent  to  18  lines 
(pairs)/mm.  Figure  15  shows  two  MTF  curves  derived  by 
scanning  bar  patterns  recorded  by  the  Dicomed  onto  2415  in 
orthogonal  directions;  the  bars  have  various  widths  and 
originate  from  the  Dicomed  operating  in  "high  resolution" 
mode.  The  modulation  caused  by  the  10  urn  aperture  employed 
in  this  microdensitometer  measurement  has  been  removed, 
along  with  the  transformation  from  square  wave  to  sine  wave 
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MTF.  Clearly,  the  Dicomed  resolution  falls  short  of  our 
intended  goal  of  reproducing,  at  high  contrast,  frequencies 
up  to  50  lin^s/mm. 


0  5  10  15  20  25 


L I NES/mm 


Figure  15:  Dicomed  MTF:  (a)  orthogonal  to  scan  direction, 
on,  (c)  average  of  a  and  b.  (b)  parallel  to  scan 
direction,  (c)  average  of  (a)  and  (b) 


No  i  se 

Fluctuations  in  density  are  evident  in  the  step  wedges 
written  by  the  Dicomed  on  the  2415  transparencies,  as  well 
as  the  SO-015  duplicates.  The  fluctuations  are  correlated 
over  a  distance  of  40  -m  (where  the  correlation  coefficient 
equals  0.1),  which  is  the  approximate  diameter  of  a  pixel 
written  by  the  Dicomed.  We  conclude  that  this  Dicomed  noise 
arises  from  fluctuations  in  the  charqe  delivered  to  the 


phosphor  when  the  Dicomed  beam  is  unblanked.  It  is  probably 
associated  with  fluctuations  in  accumulated  space  charge 
during  periods  when  the  beam  is  blanked.  This  noise 
represents  an  independent,  uniform,  additive,  uncorrelated 
contribution  to  the  noise  added  digitally.  Its  numerical 
equivalent,  in  units  of  the  digital  database,  is  given  by 
oN (Dicomed)  =  3.33.  This  number  was  determined  as  described 
in  Appendix  E. 

A  second,  less  serious  noise  was  the  apparent  outcome  of 
geometric  inaccuracy  manifest  as  dark  raster  lines  when 
viewed  under  high  magnification  (60X) ,  appearing  to  be 
spaced  at  periodically  varying  intervals  of,  typically,  100 
urn.  We  feel  that  this  form  of  "noise"  is  clearly 
discernible  as  being  a  display  defect  rather  than  being 
genuine  scene  structure.  Another  manifestation  of  geometric 
error  was  seen  intermittently  in  the  form  of  single  skewed 
lines;  one  such  example  exists  within  our  hardcopy  database 
in  scene  5200,  located  2  mm  from  the  stepwedge.  Playbacks 
whose  information  content  was  seriously  affected  by  such 
anomalies  were  repeated.  The  Dicomed  also  introduces  three 
faint  lines  perpendicular  to  the  scan  direction,  their 
spacing  being  such  as  to  divide  the  4095  x  4095  plotting 
matrix  into  four  equal  sections.  We  presume  this  to  be  an 
artifact  of  the  D/A  converter  that  controls  the  voltage 


(1)  Dicomed  geometric  accuracy  is  specified  in:  Operational 
and  Programming  Manual,  Image  Recorders,  Models  D45  and 
D47,  Dicomed  Corporation,  9700  Newton  Avenue,  South, 
Minneapolis,  MN  55431. 
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supplied  to  the  deflection  yokes  of  the  CRT. 


Tone  Reproduction 

Hardcopy  on  Kodak  2415  was  written  with  the  Dicomed  set 
in  its  "log  complement"  mode,  i.e.,  output  density  linear 
with  input  number.  The  "positive"  numerical  data  input  to 
the  Dicomed  therefore  resulted  in  a  positive  transparency. 
Tests  of  Dicomed  playback  revealed  occasional  anomalous 
behavior,  manifest  as  one  or  more  highly  non-linear  regions 
in  a  normally  near-linear  curve.  Figure  16  compares  an 
example  of  this  behavior  with  the  normal  case.  We  could 
only  speculate  as  to  the  source  of  this  disturbing  effect; 
we  surmise  that  it  is  caused  by  non-uniform  phosphor 
response  across  the  CRT  of  the  Dicomed.  The  Dicomed  has  a 
"priming  mode",  designed  to  obviate  phosphor  hysteresis; 
however,  evidence  of  questionable  post-priming  operation 
persuaded  us  to  ignore  this  feature  and  keep  the  device 
switched  on  throughout  the  film  recording  process.  Curve  C 
of  Figure  16  is  typical  of  calibrations  run  immediately 
afte-.  priming.  Curve  B  is  more  representative  of  our  data 
as  recorded  on  2415,  thouqh  some  flattening  of  the 
characteristic  at  intermediate  and  high  values,  to  a  lesser 
extent  than  that  illustrated  by  curve  C,  often  and 
unpred ic tably  occurred. 
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Figure  16:  Dicomed  tonal  transfer:  (a)  ideal,  (b) 
acceptable,  (c)  unacceptable 


Generating  Transfer  Curves  for  Playback 

Our  previously  stated  aim  for  tone  reproduction  in 
playback  was  to  maintain  the  contract  of  the  original 
Precision  Photo  positives.  Because  we  scanned  the  ground 
truth  database  in  "transmission-mode”,  there  was  a 
requirement  to  perform  a  logarithmic  transformation  on  the 
digital  database  prior  to  hardcopy  generation  on  the 
Dicomed,  whose  characteristic  is  to  yield  densities 
approximately  linear  with  input  number.  Also,  an  additional 
compensatory  transfer  curve  or  look-up  table  (LUT)  was 
needed  in  view  of  the  nonlinearity  of  the  Kodak  SO-015 
duplicating  film. 
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Design  of  LUTs  was  based  on  a  modification  relationship 
between  the  photometric  input  to  a  photographic  system  and 
the  output  of  that  system,  expressed  in  terms  of  the  final 


pho tographic 

dens i ty 

(Tupper , 

1966) . 

In  Figure  17, 

quadrants  1 

and  3 

are  known 

through 

instrument  and 

development  calibration  (quadrant  3  is  the  composite 
Dicomed/241 5/SO-Ol 5  curve)  ,  and  quadrant  4  is  the  desired 


goal  of  tone  reproduction. 


0  ,  *  D  ♦ 

out  t  out* 


It  remains  to  calculate  the  LUT  needed  in  quadrant  2  to 
provide  a  "match"  between  the  digital  database  and  the 
SO-015  display.  A  certain  output  density,  corresponding  to 
the  original  density  at  point  A,  is  traced  vertically 
downwards  to  point  B,  locating  the  equivalent  transmission 
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value  nin  in  the  digital  database.  Once  again,  starting  at 
point  A,  the  input  needed  to  produce  the  same  value  of 
output  density  on  the  playback  curve  is  located  by  tracing 
horizontally  to  point  D.  It  is  then  a  simple  matter  to  link 
points  B  and  D  to  give  a  point  C  on  the  LUT.  Additional 
points  covering  the  required  range  of  output  density  are 
then  employed  to  define  the  complete  LUT. 

The  idealized  (and  unrealizable)  case  of  a  linear 
Dicomed/contact-processing  characteristic  gives  a 
logarithmic  LUT: 


Nout  3  "out.  +  (N°ut  "  N°ut  •  )/{l0g  Nin  “  l0g  N 

min  max  min  max 


in 


min 


*  dog  Nin  -  log  N. 


in 


(10) 


min 


which,  assuming: 


0,  log  Nin  =  Jor  255, 
max 

255,  and  log:  = 

nmin 

NQut  =  (255/log  255)  log  Nin.  (11) 

In  reality,  the  composite  playback  curve  is  somewhat 
nonlinear,  as  Figure  17  shows;  actual  LUTs  were  therefore 
less  curved  than  the  above  equation  suggests.  Figure  18 
illustrates  the  three  LUTs  that  were  used,  and  identifies 


out  . 
min 


out. 


max 


g  ives : 
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the  associated  blur  levels:  (i)  2,  (ii)  0,  and  (iii)  4, 

3,1;  (i)  -  (iii)  also  represent  the  order  of  our  production 

schedule.  The  variability  between  (i)  -  (iii)  may  be 

attributed  to  two  factors:  the  fluctuation  in  Dicomed 

performance  and  USC  darkroom-processing  stability,  and  our 
increasing  understanding  of  how  the  Dicomed  photo-processing 
affected  playback  quality  and  how  to  optimize  LUT  designs  to 
compensate  for  observable  "anomalies",  such  as  gross  non- 
linearities  in  the  Dicomed  curve.  In  retrospect,  the  final 
LUT  (iii)  represents  our  best  compromise,  which,  had  the 
time  been  available,  should  have  been  used  for  blur  levels  2 
and  0  also.  Notice  that  the  dynamic  output  range  is  0  to 
243,  rather  than  0  to  255.  This  is  to  compensate  for  the 
tendency  of  the  playback  curve  to  "flatten  out"  between  243 
and  255  because  of  Dicomed  nonlinearity. 

:oe 

..  no 

.a 

100 

so 

n 

V 

o  ino  i  roc  y.n 

N. 

i  n 

Figure  18:  LUT  curves 
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SUMMARY  MEASURES  OF  QUALITY  AND  STATISTICS  OF  HARDCOPY 
DATABASE 

Blur 

An  approximate  representation  of  the  composite  MTF  for 
each  degree  of  blur,  with  no  noise  added,  is  shown  in  Figure 
19.  The  five  curves  are  derived  by  cascading  in  sequence 
three  MTF  curves:  the  PDS  MTF,  the  effecive  MTF  of  digital 

blurring,  given  by  exp  {  (-1/2)  (2Iifo)  2} ,  and,  finally,  the 
average  Dicomed  MTF. 


Figure  19:  Composite  hardcopy  MTF  at  each  blur  level. 

Solid  lines  are  ideal  (infinite  Dicomed 
resolution) 
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Table  7  compares  nominal  and  realized  hardcopy  blur 


levels . 


TABLE  7 

Nominal  and  Actual  Hardcopy  FWHM  Blur  (urn) 


Level 

Nominal 

Ha  rdcopy 

0 

20 

40 

1 

40 

52 

2 

80 

84 

3 

150 

162 

4 

320 

322 

The  disparity  is  created  by  the  relatively  poor 
performance  of  the  Dicomed  (Figure  15),  and  is  clearly  most 
serious  at  levels  0  and  1.  Hardcopy  blur  level  0  is 
approximately  the  same  as  the  nominal  value  of  40  urn  of  blur 
in  level  1. 

No  ise 

As  noted  above,  the  Dicomed  contributed  an  equivalent 
uncorrelated  additive  noise  to  our  digital  database  of  rms 
value  3.33.  Adding  this  in  quadrature  with  the  added  noise 
digitally  produced  the  actual  SNR  levels  (peak  signal  to  rms 
noise)  in  the  hardcopy  indicated  in  Table  8. 


53  - 


TABLE  8 


Nominal  and  Actual  Hardcopy  SNRs 


Level 

Nominal  SNR 

Hardcopy 

0 

200 

75 

1 

100 

60 

2 

50 

42 

3 

25 

24 

4 

12.  5 

12 

Histogram-Related  Statistics,  Tone  Reproduction 

The  hardcopy  data  shown  in  Table  9  and  the  density 
histograms  of  Figure  20  (a-j)  and  Figure  21  (a-g)  were 

simulated  on  the  International  Imaging  Systems  70/E  image 
computer  in  the  Digital  Image  Analysis  Laboratory  at  the 
University  of  Arizona.  Images  of  512  x  512  pixels,  which 
were  derived  by  subsampling  the  digital  database  at  every 
eighth  pixel,  were  stored  in  8-bit  refresh  memory  and 
hardware  LUTs  applied  in  the  IIS  70/E  pipeline  processors  in 
order  to  simulate  the  tonal  transfer  curves  and  playback 
curves  of  our  actual  hardcopy  displays.  We  have  already 
discussed  the  tonal  transfer  curves  of  Figure  18.  We 
assumed  for  our  simulation  of  hardcopy  that  the  shape  of  the 
playback  characteristic  remained  constant  as  the  maximum 
densities  varied  with  each  scene  in  the  degraded  database. 
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(By  "playback  characteristic"  we  mean  the  composite  curve  of 
the  Dicomed/2415/SO-015. )  Measurements  taken  from  hardcopy 
production  supported  the  above  assumption.  The  resultant 
density  histograms  of  Figures  20  and  21  are  essentially 
indistinguishable  from  those  of  Figures  9  and  10  with 
logarithmic  LUTs  applied.  It  should  be  noted  that  the 
abscissa  scale  of  0  255  for  Figures  9,  10,  and  13 

represents  data  values  ranging  from  0  to  255.  In  Figures  20 
and  21,  the  0  -  255  scale  is  quite  arbitrary.  It  is  a 

linear  scale  in  density,  where  0  =  and  255  =  D  .  . 

(D  .  always  lies  in  the  range  0.13  -  0.15.) 

Figure  22  compares  the  tone  reproduction  of  every  blur 
level  with  the  ideal  that  we  set  out  to  achieve,  each  curve 
being  derived  by  cascading  the  playback  curve,  LUT  (from 
Figure  18),  and  PDS  curve.  The  scale  of  the  ordinate  value 
assumes  that  we  achieved  the  correct  maximum  density.  Cases 
where  this  is  not  so  require  appropriate  correction  (see 
Appendix  D) .  Reproduction  of  contrast  for  blur  level  2  is 
obviously  the  worse  case.  A  tendency  toward  "washed-out" 
highlights  and  obscured  backgrond  details  is  the  result. 
Ironically,  some  mid-range  structure  is  enhanced,  as  in  deck 
structure  on  the  vessels  of  port  scene  51.  A  casual 
observer  would  no  doubt  notice  the  different  contrast  of 
blur  level  2  scenes.  It  remains  to  be  seen  whether  our  PI 
studies  reveal  any  disparity  in  performance  caused  by  the 
three  tonal  transfer  curves.  The  variation  in  contrast 
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among  the  remaining  blur  levels  is  not  apparent  in  our 
hardcopy,  as  perhaps  Figure  22  indicates. 


TABLE  9 

Hardcopy  Density  Histogram  Statistics 


Image  If 

Max; 

Min  Mean 

98%  D 

'.  R. 

90%  D 

>.R. 

80%  D 

.R. 

0100 

1.  85 

.15 

0.78 

1.60 

.  17 

1.60 

.  26 

1.60 

.  29 

1100 

2.  2 

.15 

1.00 

1.96 

.24 

1.71 

.33 

1.59 

.44 

4100 

2.03 

.  14 

0.  56 

1.12 

.  18 

1.05 

.  26 

0.95 

.  29 

5100 

1.88 

.15 

1.09 

1.67 

.  20 

1.47 

.34 

1.47 

.51 

6100 

1.85 

.  14 

0.  54 

1.84 

.14 

1.65 

.14 

0.82 

.16 

8100 

1.91 

.15 

1.06 

1.69 

.  28 

1.48 

.57 

1.48 

.72 

9100 

1.99 

.15 

0.94 

1.79 

.17 

1. 54 

.  33 

1.49 

.44 

1200 

1.84 

.15 

0.  86 

1.65 

.  16 

1.45 

.25 

1.  35 

.  29 

4200 

1.89 

.15 

0.77 

1.28 

.  29 

1.17 

.37 

1.07 

.42 

5200 

1. 94 

.15 

1.21 

1.73 

.  32 

1.73 

.  51 

1.52 

.63 

0120 

1.9 

.15 

0.84 

1.82 

.17 

1.80 

.  22 

1.80 

.  25 

4120 

1.79 

.  14 

0.55 

1.55 

.19 

1.43 

.  22 

1.  25 

.  25 

9120 

1.63 

.15 

0.86 

1.53 

.16 

1.53 

.55 

1.45 

.  26 

0140 

1.  85 

.15 

0.  80 

1.65 

.  17 

1.65 

.  26 

1.65 

.29 

4140 

1.41 

.15 

0.62 

1.19 

.  26 

1.  12 

.  30 

1.02 

.  34 

9140 

1.75 

.  16 

0.96 

1.65 

.  18 

1.63 

.37 

1.45 

.45 

« 
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MAX  -  25*1,  MIN  *  0,  MEAN  *  1A0 

(b)  0140 


MAX  =  255,  MIN  =  36,  MEAN  - 

(c)  4120 


(d)  4140 


U)  9120 


(f)  9140 


Figure  21:  Representative  hardcopy  histograms  (Blurs  2,  4: 

Noise  0).  Abscissa  scale:  0  =  D^,  255=  Dmin  (Dmin“ 
0.13  -  0.15.) 

Values  of  D^; 

(a)  1.9,  (b)  1.85,  (c)  2.1,  (d)  1.79,  (e)  1.7,  and  (f)  1.83. 
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(ii)  BLUR  0 
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ORIGINAL  FILM  DENSITY 

Figure  22:  Hardcopy  tone  reproduction. 


Table  10  given  the  measured  maximum  density  of  all  250 
scenes  from  the  selected  dark  regions  described  in  Table  11. 


TABLE  10 


Hardccpy  Maximum  Scene  Densities 


Noise  Level 

0 

1  2 

3 

4 

010 

1.58 

1.85 

1.58 

1.27 

1 .  23 

Oil 

1.74 

1.73 

1.55 

1 . 56 

1.45 

012 

1.78 

1.3  2 

1.52 

1.44 

1 . 26 

013 

1.71 

1.78 

1 . 54 

1.63 

1.35 

014 

1.74 

1.72 

1.55 

1.64 

1.45 

110 

1.50 

1.70 

1.89 

1.73 

1.  57 

111 

1.50 

1. 53 

1.88 

1 . 56 

1.47 

112 

1.98 

2.07 

1.98 

2.  21 

1.85 

113 

1.41 

1 . 48 

1.49 

1.69 

1.69 

114 

1.50 

1.53 

1.53 

1.77 

1 .  50 

410 

0.97 

1.23 

0.89 

1.00 

1.12 

41.1 

1.24 

1.30 

1.31 

1.40 

n/a 

412 

1.59 

1.41 

1.25 

1.55 

1.26 

413 

1.23 

1.24 

1.28 

1.30 

1.20 

414 

1.15 

1.18 

1.00 

1.23 

1.10 

510 

1.82 

1.81 

1.73 

1.60 

1.15 

511 

1.32 

1 .  43 

1.48 

1.52 

1.38 

512 

2.22 

2.04 

2.31 

2.10 

2.09 

513 

1.52 

1.57 

1.53 

1.44 

1.31 

514 

1.74 

1.71 

1.59 

1.  59 

1.51 

510 

0.59 

0.74 

0.88 

0.90 

1.10 

511 

1.44 

1 . 44 

1.4  8 

1.53 

1.44 

512 

1.54 

1 . 88 

1.51 

1.15 

1 .  50 

513 

1.08 

1.11 

1.25 

1.25 

1.38 

514 

1.01 

0.98 

0.98 

1.02 

1.16 

310 

1.88 

1.70 

1.82 

1.55 

1.  54 

811 

1.54 

1.52 

1.50 

1.37 

1.27 

812 

3.0  0 

2.88 

2.  37 

1.45 

2.05 

813 

1.  54 

1.44 

1.35 

1 .  18 

1.18 

814 

1.70 

1.85 

1.50 

1.50 

1.50 

910 

1. 80 

1.72 

2.02 

2.03 

1.51 

911 

1.82 

1.55 

1.53 

1.56 

1.45 

912 

1.38 

1.52 

1.62 

1.49 

1.46 

9l  3 

1.78 

1.88 

1.70 

1 . 62 

1.33 

914 

1.80 

1.77 

1.71 

1.59 

1.57 
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120 

1.53 

1.49 

121 

1.44 

1.47 

122 

2.56 

2.60 

123 

1.21 

1.35 

124 

1.36 

1.39 

420 

1.23 

1.33 

421 

1.42 

1.41 

422 

1.15 

1.74 

423 

1.23 

1.35 

424 

1.30 

1.35 

520 

1.54 

1.71 

521 

1.59 

1.61 

522 

0.91 

0.90 

523 

1.67 

1.63 

524 

1.42 

1.54 

1.31  1.36  1.16 

1.45  1.42  1.42 

2.37  2.32  1.88 

1.38  1.51  1.32 

1.33  1.47  1.49 

1.39  1.48  1.44 

1.51  1.58  1.61 

2.04  2.42  1.85 

1.47  1.32  1.45 

1.35  1.40  1.30 

1.78  1.50  1.17 

1.55  1.56  1.44 

1.00  0.99  0.69 

1.43  1.42  1.31 

1.56  1.60  1.44 


TABLE  11 

Areas  of  Maximum  Density,  by  Scene 

Scene  Description  of  Area  Where  Density  Was  Recorded 

01  Shadow  of  largest  building 

11  Shadow  of  largest  building 

41  Sea  near  vicinity  of  ship  bows 

51  Most  open  area  of  sea  in  corner  of  scene 

61  Shadow  of  building 

81  Dark  area  adjacent  to  third  darkest  step  in 

stepwedge 

91  Largest  area  of  shadow  (shadow  of  building) 

12  Shadow  of  building 

42  Open  sea  in  empty  corner 

52  Open  sea  in  corner  opposite  stepwedge 
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Superimposed  on  the  variability  due  to  photographic 

processing  is  a  number  of  trends  associated  with  the  amount 
of  noise  present.  For  example,  if  the  darkest  scene  area 
has  signal  level  0,  then  additive  noise  will  result  in  a 
skewed  Gaussian  distribution: 

p(Z)  =  l/(d5Tir)  exp  (-(£2)/(2 a2)).  (12) 

The  expected,  or  mean  level  £  will  be  that  measured  by  the 
densitometer,  since  it  has  an  integrating  effect,  i.e., 

-  r  i  .255 

£=[l/(af2pj  j  £  •  exp(- (£2)/(202)  )  d£  ,  (13) 

o-' 

which  approximates 

£  *  [l/(  afTrr)J  a2  =  o.4a.  (14) 

Our  four  levels  of  noise  give  expected  values  1,  2,  4, 

and  8,  which  typically  result  in  densities  of  1.68,  1.60, 

1.49,  and  1.14  in  hardcopy.  This  reduction  of  background 
density  with  increasing  noise  is  evident  in  scenes  01,  51, 

81,  and  91.  Interestingly,  the  airfield  scene  61  exhibits 
the  reverse  trend.  This  is  due  to  the  high  mean  brightness 
(low  mean  density)  apparent  in  Table  9  and  the  broad 
histogram  of  Figure  20  (e)  . 


VI.  DISCUSSION 


We  have  reported  our  attempts  to  develop  an  imagery 
database  from  high  quality  aerial  photography  relevant  to 
the  image  interpretation  task.  Specifically,  we  have 
created  two  versions  of  a  database  of  250  images,  consisting 
of  10  scenes  each  represented  at  5  levels  of  noise. 

The  digital  database  is  in  8-bit  format  on  9-track  tape 
with  4095  bytes  (pixels) /record  and  4096  records/image  file. 
We  have  achieved,  to  high  accuracy,  the  specifications  for 
blur  and  additive  noise.  The  effect  of  signal-dependent 
noise  caused  by  the  grain  in  the  original  photographic 
transparencies  needs  to  be  studied.  At  present,  all  data 
derived  from  digitization  of  the  tr ansparencies  are  regarded 
as  "signal",  or  "ground-truth". 

The  hardcopy  database,  on  Kodak  2415  film,  suffers  from  a 
number  of  drawbacks  not  seen  in  the  digital  version. 
Problems  were  encountered  in  the  control  of  film-playback  on 
the  Dicomed  and  in  darkroom  repeatability.  Nevertheless, 
these  variations  are  well  quantified  in  this  report  and  the 
hardcopy  database  is  useful  and  representative  of  our 
operational  image  quality.  In  spite  of  the  limitations  of 
the  hardcopy  imagery,  very  good  PI  data  were  obtained, 
supporting  the  notion  that  errors  in  the  database  were  small 


compared  to  the  main  effects  of  the  "ariables  of  interest. 
These  results  are  described  in  other  reports  in  this  series. 

The  digital  database  is  clearly  better  controlled  than 
the  hardcopy  database.  We  feel  that  since  most  of  our 
program  will  use  the  digital  database  in  subsequent 
experiments,  then  the  reduced  consistency  of  hardcopy  is  not 
particularly  critical. 
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Appendix  A 

NATO  IMAGERY  INTERPRETABI LITY  RATING  SCALE 

The  following  is  an  extract  from  the  Image 
Interpretabil ity  Scale  published  by  the  NATO  Air 
Standardization  Coordinating  Committee  (ASCC-AIR  STD  101/11, 
10  July  1978) .  It  is  very  similar  to  the  heavily  employed 
National  Image  Interpretabil ity  Rating  Scale  (NIIRS). 
Rating  categories  0  through  3  are  not  presented,  since  they 
represent  ground  resolved  distances  poorer  than  2.4  m.  Such 
low  resolution,  according  to  photo  interpreters  employed  in 
our  studies,  is  rarely  encountered  in  aerial  photographs  of 
tactical  utility. 

The  terms  "detection",  "recognition",  "identification", 
and  "technical  analysis"  are  also  defined  herein.  Each  term 
represent  an  increasingly  demanding  task  requiring 
increasing  resolution.  Hence,  "technical  analysis",  for 
example,  is  prevalent  in  the  descriptions  of  rating 
categories  8  and  9. 

DEFINITIONS 

G round  Resolved  Pis tances ;  Ground  resolved  distance  is 
the  minimum  test  target  element  resolved  on  the  ground. 
With  a  system  that  produces  a  Ground  Resolved  Distance  (GRD) 
of  0.3  m,  the  smallest  bar  of  the  test  target  that  is 
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distinguishable  has  a  physical  width  of  0.15  m.  (A  tri-bar 
test  target  was  used  to  determine  GRD  and  subsequently  to 
calibrate  the  Imagery  Interpretabil ity  Rating  Scale). 

Detection:  In  imagery  interpretation,  the  discovering  of 

the  existence  of  an  object  but  without  recognition  of  the 
object. 

Recogni tion :  The  determination  by  any  means  of  the 

friendly  or  enemy  character  of  the  individuality  of  another, 
or  of  objects  such  as  aircraft,  ships,  tanks,  or  of  a 
phenomenon  such  as  communications  or  electronics  patterns. 

Identification:  In  imagery  interpretation,  the 

discrimination  between  objects  within  a  particular  type  or 
class . 

Technical  Analysis :  The  ability  to  describe  precisely  a 

feature,  object,  or  component  imaged  on  film. 

Rating  Category  4. 

*  Ground  Resolved  Distance,  1.2  to  2.5  m.  (3.94  to  8.2 
ft.)  (47.  25  to  98  in.) 

*  Detect  rockets  and  artillery. 

*  Recognize  troop  units. 

*  Recognize  aircraft  (such  as  FAGOT/MIDGET  when  singly 
deployed)  . 

*  Recognize  missile  sites  (SSM/SAM) .  Distinguish  between 
missile  types  by  the  presence  and  relative  position  of 
wings  and  control  fins. 

*  Recognize  nuclear  weapons  components. 


*  Recognize  land  minefields. 

*  Identify  ports  and  harbors. 

*  Identify  railroad  yards  and  shops. 

*  Identify  trucks  at  ground  force  installations  as  cargo, 
flatbed,  or  van. 

*  Identify  a  KRESTA  by  the  helicopter  platform  flush  with 
the  fantail,  a  KRESTA  II  by  the  raised  helicopter 
platform  (one  deck  level  above  fantail  and  flush  with 
the  main  deck) . 

Rating  Category  5. 

*  Ground  Resolved  Distance,  0.75  to  1.2  m.  (2.46  to  3.94 
ft.)  (29.53  to  47.  25  in.) 

*  Detect  the  presence  of  call  letters  or  numbers  and 
alphabetical  country  designator  on  the  wings  of  large 
commercial  or  cargo  aircraft  (where  alphanumerics  are 
three  feet  high  or  greater)  . 

*  Recognize  command  and  control  headquarters. 

*  Identify  a  singly  deployed  tank  at  a  ground  forces 
installation  as  light  or  medium/heavy. 

*  Perform  Technical  Analysis  on  airfield  facilities. 

*  Perform  Technical  Analysis  on  urban  areas. 

*  Perform  Technical  Analysis  on  terrain. 


*'  * 
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Rating  Category  6. 

*  Ground  Resolved  Distance,  40  to  75  cm.  (1.31  to  2.46 
ft.)  (15.75  to  29.5  in.) 

*  Recognize  radio/radar  equipment. 

*  Recognize  supply  dumps  (POL/Ordnance) . 

*  Recognize  rockets  and  artillery. 

*  Identify  bridges. 

*  Identify  troop  units. 

*  Identify  coast  and  landing  beaches. 

*  Identify  a  FAGOT  or  MIDGET  by  canopy  configuration  when 
singly  deployed. 

*  Identify  the  following  ground  force  equipment:  T-54/55 
tank,  BTR-50  armored  personnel  carrier,  57  mm  AA  gun. 

*  Identify  by  type,  RBU  installation  (e.g.,  2500  series), 
torpedo  tubes  (e.g.,  21  in./53.34  cm),  and  surface-to- 
air  missile  launchers  on  a  KANIN  DDG ,  DRIVAC  DDGSP,  or 
KRESTA  II. 

*  Identify  a  ROMEO-class  submarine  by  the  presence  of  the 
cowling  for  the  snorkel  induction  and  the  snorkel 
exhaust. 

*  Identify  a  WHISKEY-class  submarine  by  the  absence  of 
the  cowling  and  exhaust. 

Rating  Category  7. 

*  Ground  Resolved  Distance,  20  to  40  cm.  (0.66  to  1.31 
ft.)  (7.87  to  15.75  in.). 

*  Identify  radar  equipment. 
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*  Identify  major  electronics  by  type  on  a  KILDEN  DDGS  or 
KASHIN  DLG. 

*  Identify  command  and  control  headquarters. 

*  Identify  nuclear  weapons  components. 

*  Identify  land  minefields. 

*  Identify  the  general  configuration  of  an  SSNB/SSGN 
submarine  sail,  to  include  relative  placement  of  bridge 
periscope (s)  and  main  electronics/navigation  and 
equipment. 

*  Perform  Technical  Analysis  on  ports  and  harbors. 

*  Perform  Technical  Analysis  on  railroad  yards  and  ships. 

*  Perform  Technical  Analysis  on  roads. 

Rating  Category  8. 

*  Ground  Resolved  Distance,  10  to  20  cm.  (0.33  to  0.66 
ft.)  (3.94  to  7.87  in.) 

*  Identify  supply  dumps  (POL/Ordnance) . 

*  Identify  rockets  and  artillery. 

*  Identify  aircraft. 

*  Identify  missile  sites  (SSM/SAM) . 

*  Identify  surface  ships. 

*  Identify  vehicles. 

*  Identify  surfaced  submarines  (including  components  such 
as  ECHO  II  SSGN  sail  missile  launcher  elevator  guide 
and  major  electronics/navigation  equipment  by  type) . 

*  Identify,  on  a  KRESTA  II,  the  configuration  of  major 
components  of  larger  electronics  equipment  and  smaller 
electronics  Dy  type. 
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*  Identify  limbs  (arms,  legs)  on  an  individual. 

*  Perform  Technical  Analysis  on  bridges. 

*  Perform  Technical  Analysis  on  troop  units. 

*  Perform  Technical  Analysis  on  coast  and  landing 
beaches. 

Rating  Category  9. 

*  Ground  Resolved  Distance,  less  than  10  cm.  (<  0.33  ft.) 
(<  3.94  in.) 

*  Identify  in  detail  the  configuration  of  a  D-30  howitzer 
muzzle  brake. 

*  Identify  in  detail  on  a  KILDEN  DDGS  the  configuration 
of  torpedo  tubes  and  AA  gun  mountings  (including  gun 
details) . 

*  Identify  in  detail  the  configuration  of  an  ECHO  II  SSGN 
sail  including  detailed  configuration  of  electronics 
communications  equipment  and  navigation  equipment. 

*  Perform  Technical  Analysis  on  radio/radar  equipment. 

*  Perform  Technical  Analysis  on  supply  dumps 
(POL/Ordnance) . 

*  Perform  Technical  Analysis  on  rockets  and  artillery. 

*  Perform  Technical  Analysis  on  missile  sites. 

*  Perform  Technical  Analysis  on  nuclear  weapons 
components . 
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Appendix  B 

ILLUSTRATION  OF  DATABASE 

For  purposes  of  database  demonstration,  we  use  an  aerial 
scene  of  the  U.S.  Army  Tank  and  Automotive  Command  (TACOM) , 
Warren,  Michigan.  (This  scene  was  not  included  in  the  10 
photographs  of  the  database.)  The  original  negative  was 
exposed  at  an  altitude  of  1820  ft.  with  a  lens  of  6  in. 
focal  length,  yielding  a  scale  of  1:3640.  Digitization  at 
20  ym  intervals  is  therefore  equivalent  to  a  "ground  truth" 
resolution  of  approximately  2.9  in.  Figure  23  shows  this 
ground  truth  image,  with  no  additional  blur  or  noise,  in  a 
form  illustrating  only  every  eighth  pixel  of  every  eighth 
line  out  of  the  total  of  4096  x  4096  pixels. 

As  described  in  the  text,  we  generated  the  following 
equivalent  GRD  values,  by  means  of  digital  blurring: 
"ground  truth",  15,  30,  61,  and  122  cm.  Figure  24 
illustrates  a  matrix  of  nine  degraded  images,  covering  an 
area  of  512  x  512  pixels  at  full  resolution,  selected  from 
the  global  scene  of  Figure  23.  We  show  only  three  levels  of 
blur  (ground  truth,  30,  and  122  cm)  and  three  levels  of 
noise  (nominal  SNRs  of  200,  50,  and  12.5). 
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Appendix  C 


SCAN  CHARACTERISTICS  AND  SETUP  PROCEDURE  FOR  THE 
PDS  MICRODENSITOMETER 


The  original  digitized  data  were  prepared  by  scanning  the 
Precision  Photo  positive  transparencies.  The  scanning  was 
performed  on  a  PDS  Model  1010A  microdensitometer.  The 
specifications,  as  supplied  by  the  manufacturer,  are  listed 
in  Table  12.  Figure  25  shows  a  rough  sketch  of  the  optical 
system  layout  for  this  instrument. 

Data  were  generated  by  scanning  in  transmission  mode, 

L.e.,  N  =  KT  where  K  is  some  constant.  The 
out 

microdensitometer  is  capable  of  generating  12  bits  of  data 
per  pixel.  For  convenience,  the  scans  were  set  up  so  that 
the  maximum  transmission  reading  for  a  particular  scene 
would  yield  a  digital  value  of  2000.  There  are  three 
adjustments  to  be  made  when  setting  up  the  scans  to  insure 
that  the  data  will  lie  within  this  range  of  0  to  2000. 
First,  the  region  of  the  film  with  the  highest  transmission 
is  placed  in  the  beam.  Then,  various  neutral  density 
filters  are  inserted  in  the  beam  until  a  specified 
photomultiplier  tube  (PMT)  current  has  been  reached.  If 
this  PMT  current  cannot  be  obtained  by  changing  neutral 
density  filters  alone,  the  filter  which  comes  closest  to 
yielding  this  specified  value  is  left  in  the  beam.  At  this 
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TABLE  12 


Specifications  of  PDS  1010A 


PDS 

Microdensitometer : 

Specifications: 

Stage  Motion 

X  and  Y  travel 

Positioning  Accuracy 

Straightness  of  Travel 

254  mm  max 

+1  pm 

+5  Pm 

Illumination 

System 

Microscope 

Microscope 

Apertures 

Objective 

Eyepieces 

10X,  0.25  NA 

5X,  10X,  15X,  20X 

8  Selectable  Apertures 

Scanning 

System 

Microscope 

Microscope 

Apertures 

De  tec tor 

Objective 

Eyepieces 

10X,  0. 25  NA 

5X,  10X,  15X,  20X 

8  Selectable  Apertures 

Sll  PMT,  25  mm  diameter 

stage,  the  voltage  on  the  PMT  dynodes  is  varied  until  the 
desired  output  current  is  obtained.  Now,  two  further 
adjustments  need  to  be  made.  The  voltage  generated  by  the 
PMT  current  is  amplified  before  being  input  to  the  A/D 
converter.  The  gain  of  this  amplifier  is  adjusted  so  that 
the  signal  input  yields  a  digital  value  of  2000.  Also,  the 
zero  level  of  this  differential  amplifier  must  be  adjusted 
so  that  the  dark  current  of  the  PMT  gives  an  output  value  of 
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Figure  25:  Preslit  light  path 


zero.  This  is  done  by  blocking  the  beam  with  an  opaque 
filter  and  adjusting  the  offset  voltage  of  the  amplifier. 
The  gain  and  offset  adjustments  for  the  amplfier  must  be 
iterated  until  both  conditions,  MAX  =  200,  MIN  =  0,  are  met. 

Figure  3  of  this  report  illustrates  the  relationship 
between  the  number  output  by  the  A/D  converter  and  the 
diffuse  density  of  the  area  being  scanned.  Although  the 
exact  intercept  along  the  log  N  axis  is  weakly  dependent 
upon  the  gain  of  the  amplifier,  the  only  major  differences 
among  the  curves  are  their  lateral  displacement  and  their 
nonlinear  tails  for  small  values  of  log  N.  Over  the  8-bit 
range  of  interest  here,  the  curves  are  well  approximated  by: 

D  -  i,Dmin  '  <W/l09  W>  N  +  Dmax  <cl> 
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The  shape  of  the  nonlinear  tail  of  the  log(N)  versus 

Ddiffuse  curve  arises  as  follows.  The  linear  portion  of  the 
curve  suggests  that  the  scanning  geometry  of  the 
microdensitometer  is  such  that  T  ^  =  (t,.,,  (The 

specific  slope  depends  on  the  film  and  the  numerical 
apertures  of  the  densitometer.)  Now  the  current  from  the 

PMT  consists  of  a  contribution  due  to  the  signal,  plus  a 

1  40 

contribution  due  to  dark  current,  i.e.,  IpN(T  =  ^(T^jr  *  + 

0DC)  where  0  is  the  dark  current  term.  This  current  is 
converted  to  a  voltage  and  amplified  by  a  linear 
differential  amplifier  before  being  fed  to  an  A/D  converter. 
Thus, 


Vout  =  K2{  (K  1  (Td  i  f  +  °DC)  +  Voff*’ 


(C2) 


Here  K?  is  a  constant  which  represents  the  gain  setting  of 
the  amplifier,  while  is  the  reference  voltage  level 
which  is  adjusted  so  that  when  only  dark  current  is  present. 


V 


out 


0, 1  .e  . ,  VQff  K1  0DC. 


Since  the  number  recorded  by  the  microdensitometer  is 
proportional  to  the  output  voltage,  we  can  see  that 


Nout  -  <K' 


DC  +  K2Voff' 


(C3) 


If  the  offset  is  accurately  set. 


log  Nout  =  log  K'Tdif 


1.  40 


( C4 ) 


If  there  is  some  error  in  the  offset  setting,  then 
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log  Nout  =  log(KiT 


+  £  )  • 


(CS) 


1.40 

Equation  (C4)  implies  that  the  log  N  versus  is  a 
straight  line  of  slope  -1.40.  Equation  (C5)  predicts  the 
same  behavior  for  large  values  of  Td-ff  provided  is  small. 
However,  for  small  Equation  (C5)  predicts  a  tail  whose 
exact  shape  depends  upon  the  sign  and  magnitude  of  e. 

The  leftmost  curve  of  Figure  3  illustrates  the  effect  of 
a  (typical)  small,  negative  offset  error.  The  other  two 
curves  exhibit  the  (intentionally  produced)  effects  of 
larger  positive  offset  errors. 
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Appendix  D 


TECHNICAL  SPECIFICATIONS  MODEL  D47  DIGITAL  COLOR 

IMAGE  RECORDER 
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t*)  - ,s 7/ie  Model  D47  Digital  Color  Image  Recorder  converts 
/■:  'rqf&P  data  to  high-resolution  color  or  black  and  white  film.  ,.l'-i-* 

;,/f •'■  Each  da /a  word  (exposure  code)  transmitted  to  the  D47  producMsWa^ret 
fgl\4i  picture  element  (pixel).  Each  pixel  is  constructeeUeLoae  or  mor&WtbQ&Ui 
tT*  ■’'..‘‘■spots  generated  by  focusing  a  CRT  (cathode-ray  tub^Z&Slf^Vhtd^tbBiJjJnl^^SKW 
a-. plane.  Exposure  energy  is  either  linearly  or  logarithmically  related  to  thefjfflffi 
'vjmif*’,  date  word  value.  Jam  —  '  ’ ‘'v 
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■'}■  posure  through  the  various  filters  produces  a  particular  color  depeAdlHg^ts^J 
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The  047  can  operate  on-line  from  a  digital  computer  or  off-line  from" 
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Bulletin  12G123— Illustrated  brochure  describes  basic  recorder  oper-^x-zXtyC 
ation,  overall  control  flexibility ,  and  plotting  ver-,  tility  of  film  recorders./^'^f 
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'  •<  Eacrt  data  word  (exposure  code)  transmitted  to  the  D47  producdl&Wsbrei 
picture  element  (pixel).  Each  pixel  is  constructe&eLone  or  moiewtbQptif 
Wk  h-k-  spots  generated  by  focusing  a  CRT  (cathode-ray  tubcj^SfTTOntchUi^jJr^ 
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f  'v*  Color  recording  is  achieved  by  the  sequential 

three  (blue,  green  &  red)  spectral  filters  in  the  CRT  bea/frpefnr^l^mS 
1  < ;  posure  through  the  various  filters  produces  a  particular  color  depend, 
v  on  the  relative  exposure  energy  levels  tor  each  liltt  •- 

'  77te  D47  can  operate  on-line  from  a  digital  comput  er  or  off-line  frorn  'm&g 

inetic  tape.  Operation  under  full  program  control  is  achieved  v 
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-  Bulletin  1 2G1 23— Illustrated  brochure  describes  basic  recorder  oper-]i'» 
‘  ation,  overall  control  flexibility,  and  plotting  vers  utility  of  film  recorders..^ 

2'V  Bulletin  1 2DQ34  —  Technical, Specifications  lor.  1 1odeiD46  Black  /White';# 
4‘  Film  Recorder. 
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GEOMETRIC  PERFORMANCE 


“The  geometric  specifications  shown  below  cover  three  different  film  formats  as  well  as  two  different 4 
CRT'Configurations;  high  resolution  CRT  (standard)  and  uitra-high  resolution  'CRT  (optional).  wjjj 


EFFECTIVE  RECORDING  AREA: 


PIXEL  CONSTRUCTION:  Each  picture  element  (pixel)  is  a 
constructed  by  one  or  more  CRT  spot  exposures  depending  J 
on  the' plotting  resolution  selection  (see  Tigure)^^,^s^.^| 

High  Resolution:  1  spot  per  pixel 

Medium  Resolution:  4  spots  per  pixel 

Low  Resolution.1 16  spots  per  pixeL-jjffiBS^PI^f^^pl 


Kp.UUUt. 

i(P*se,.'<  ,v  lows 
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PLOTTING  MATRIX:  Each  D47  offers  push-butfon  or  pro-'^i 
gram  selection  of  three  different  plotting  matrices  yielding 
three  different  resolutions.  •'  '  ; 


Lowp^oiUM 

Kf v  "  *  PIXEL  CONSTRUCTION^^! 
Kf  FOR  LOW  MEDIUM  AND  HIGH  ™ 
hf  PLOTTING  RESOLUTION.  SPOT.SIZE 
IS. .CONSTANT.  <$$$&£& 


*■  ^  ’  STANDARD  CRTTOPTIONAL  CRT 

S-  Low  Resolution  512x512  pixels  £1024x1024 pixels 
L  Medium  Resolution.  1024x1024  pixels *2048x2048 pixels''^Jpj£wT ]» 
JLHigh  Resolutions  ^  2048x2048  pixels  .4096x4096  pixels  ;4-"T 


Medium  Resolution .  1024x1024  pixels  2048x2048  pixels'® )» . 
JLHigh  Resolutions.^  2048x2048  pixels  .4096x4096  plxels^^^^^" 

RESOLUTION:  Standard  High  Resolution  CRT:  2000  re- 
solvable  lines.  Optional  Ultra-high  Resolution  CRT:. 3000 
resolvable  lines.  '  -  ••.  ? 

NOTE:^Re solution  measured  on  Plus  X  Pan  .4^^film^j^^gE^j^^' 

GEOMETRIC  ACCURACY:  "Percent'1  specifications  are  ref^«W||jSjiw& 
erenced  to  the  major  axis  (Hc  or  Vc )  of  the  recording  ma-ra(|^g|jBc'. 
trix.  Refer  to  the  figure  below  for  definitions.pf/ ecojdtng  raMXb'. 

2;' matrix  measurements. 


.IV.1  V  -SiJ 


OKTf  ft 

•  VyJ  4).: 


■mm;  r:>, 


?.=r  r-.Hc-5^ 


Wifi.  Matrix  Orthogonality  (0);  ±0.25%  max.ffigjry 

£  *•*'  \Mt+M}  /  v 


mm 


X  Matrix  Rectangularity  (R):  ±0.3%  max.^^^J 
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*»»•  HIUIMA  MOVtUM^U'UI  ItJ  »/■  - V. V  II  IUA.  MBR] 

«■»  .  ^  /  (j  ..  v  i  *y»rjy?^MS 

•>:« .  ;  fl=±r-gi^-)x  ioo-"^^m| 

\  Hc+Vc  /  .  ')  '••It* 

Jj?»  Line  Curvature  (P);  {pin  cushion  distortion):  ±0.15%  maxs 

;>  ptL  (  d' \  100  or  P-  +  (£-)x  100  vf  M 
If;;  ■  iff  \Hc /  \Vc  /  ■  I  “7 


h2 


ft*** 


•sri"'. 


>1:..  *4 

:  A1.,  sj 


x  100  v 


fy.  '  d  -  Maximum  deviation  from  best  fit  straight  line. . 
V.  Spaclal  Repeatability  (S):  ±0.05%  max.  •  •*  ££  y 

i  -  1  /h\  _  /rt\ 


?•*(£) 


x  100  or  s  -  ±  100  ■  ^ 


d  =  Maximum  .displacement  of  any  given  150101 tor^ 
^successive  images  repeated  during  a  20-minute 
'^interval  .."  r  s-J 


BMP  r  PLOTTING  MATRIX  T  J 

pg^,  MEASUREMENT  DEFINITIONS/^ 

gK8|^  ARRAY  GEOMETRICS 

Jjjjp  H,  =  Length  of  top  horizontal  line  j*  .JyT'  ^ 
Ky  Hc  *  Length  of  horizontal  maiot  axis  ’  Ai 
3?  H2  -  Length  ot  bottom  horizontal  line  T.:' 

Bv  V,  =  Length  of  left  vertical  line  ''•.'tVl :'•'  'H 

k'|,  Vc  =  Length  of  vertical  major  axis  V:-  -rAf.'.oJ 

TWt  V2  -  Length  of  right  vertical  line  '  % 

M)  *  Length  of  0,0  to  max.,  max.  matrix  diagonal 
|fP„  M2  =  Length  of  0,  max.  to  max.,  0  matrix  diagonal m 


Direction  of  scan  line 


^0,n,  Spacing  Linearity  (L):  ±  0  35%  Max.  '  "  7  *  ™ 
f'  u  =1  Snm^x-S^mm  ^,0DU=  ^  Swnax-Svmin  x  J(K) 
‘  ■'*  "  2Hd  2Ve  •  ',f 

./iih  Smax  and  Smm  are  the  maximum  and  minimum  re-^, 

/  K  sultant  line  segments  constructed  from  an  internally^ 


'  V-  tj 


m 


RECOMMENDED  FILM  TYPES 


L\  r’-’  7jm  ^  «.  r-  .  ;  £•  jjtiJKjjif  ttiacx.  wmie  iiansoaienr 

t  'ff  Th  4±|-£— X'  100  ^  f'PW  Bl.ck/Whita  Trinspirew 
X?  ■  V L  Hl  +Hz  J  .  Colo.  Print 

N  >  •>  VCrlJr^1  1~*  y  y  I  '  Blatlt/Wh.le  Print 

t\r  a?  ‘  VattTv  =  ±  - - — x  100  -3ueUiUiU4<ML  ti  B'atkWMe  Print »  >. 

I'-’  '  '■■QfojJ?-'  L'*  It  +44  <*  j  >  .  y-m,n 


USE 

{Positive  Color  Transparency 

(Positive  Color  Transparency 
f  Vesaf  ■ 

‘  Black /While  Tiansoaienry 
Black/Whits  Transparency 


v  FORMAT  *r'  FILM  TYPE 

4  x  5"  sheet  •’■x’  Kodak  fktachiome  61 1 

70ram,  105mm  roll  -i-.  Kodak  Eklachiome  MS 
' * Ae.opr3phic  2448 
4  <  5M  sheet  K°dak  Super-XX  4142 


4"  i  5"  sheet  Kodak  Super*XX  4 

4"  x  5"  sheet  Kodak  Plus  X  Pah 


4  x  5"  sheet  . ' '  jj  Polaroid  Type  58 
4"  x  5"  sheet  ’y^.  Polaroid  Tyoe  52 
4"  x  5”  sheet  T-  po!i!0'!*  Tvp;  ct 


Copy  avoilable  to  DTIC  does  not 
permit  fully  legible  reproduction 
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PHOTOMETRIC  PERFORMANCE  ■  W 


The  D47  is  capable  of  converting  an  exposure  code  into  any  one  of  64  exposure  levels  (6-bit  codd)  on 
256  exposure  levels  (8-bit  code).  The  constant  intensity  CRT  beam  is  time-modulated  to  produce, ex 
posures  either  linearly  or  logarithmically  related  to  the  exposure  code.nv 


EXPOSURE  RANGE:  Depends  on  film  type;  range  is  2  OD  EXPOSURE  CONTROL:  The  CRT  beam  intensity  may  be,1 
(diffuse  Density  units)  for  Ektachrome  6115  film  and  1.8D-  varied  to  compensate  for  different  film  speeds.' Note: 
fcrPlusX  Pan  4147  film.  ’  The  exposure  level  may  be  periodically  checked  and  re-1 

cYonciiDc  iiMieoaMrrv*  calibrated  by  use  of  the  digital  panel  meter  on  the  TEST, 


uniformity  Module  (HUM).  ^ 

STANDARD  UNIFORMITY:  ±0.0750  max.'  % 
UNIFORMITY  WITH  HUM:±0.035D  max,,  '"  ^ 


RECORDING  SPEED 


.  •  <•  ;Vv**  \  .  •  « 

The  total  time  required  to  record  an  image  is  a  function  of  (1)  plotting  matrix  size  (number  of  CRT  spots), 

.  (2)  the  average  exposure  level  required,  and  (3)  the' number,  of  filters  selected.'-The  time  required  td 
'  record,  any  given  pixel  is  a  function  of  (1)  number  of  spots  per' pixel  and  (2)  the  exposure  level  required^# 

■  ?  •  v  <t  •«*  V  •  ••  /ixarNt-M 

SPOT  EXPOSURE  TIME:  7^  sec.  average  code;  30,'  sec.’ IMAGE  GENERATION  TIME:  The  following  limes  are* 
max.  NOTE:  times  are  for  recording  in  leg  mode JBwmS&tt. those  fequired  for  recording  the  full  matrix  inMog  mode* 


I/O  TIME:  2/r  sec.  max. 


ADJACENT  SPOT  POSITIONING 


RANDOM  POSITIONING  TIME:  1.0  msec,  per  axis  optional  CRT:’. 5.5  minutes  in  black/white;  16.5  minutes  in1 

FLYBACK  TIME:  1 .0  msec.  ' v .'  . 


COLOR  SYSTEM 


1 


Full  color  recording  is  achieved  by  sequentially  exposing  the  film  through  a  series  of  compensated ^ 
spectral  filters.  The  entire  image  is  first  exposed  through  one  filter  and  the  next  filter  is  rotated  into) 
position  (program  or  push-button  commanded)  tor  the  next  exposure.  A  blank  position  is  provided  tor  , 
black/white  recording.  ‘ 


BLUE  FILTER: 


Wraiien  #4/ 


■  * :  i  ’  r 


GREEN  FILTER 


Wratten  458B 


RED  FILTER 


Wratten  425 


CRT  "’HOSF’HOR: 


•  /><  '■  V-'bL'h.  , 


OPERATING  &  TEST  PANELS 


l 


The  recorder  includes  an  operations  control  panel  and  a  test  (and  calibration)  pane',  items  on  the  test 
panel  are  indicated  by*.  1  ”  '  ■  •/•/*  * 

POWER:  Key  operated  OFF,  ON  &  LOCK-ON  positions  with  FILTER:  Push-button  selection  of  NEUTRAL,*  RED.  BLUf 
power  indicator.  .  or  GREEN  filter. 

INPUT:  6-BIT  or  8-BIT  exposure  code  push-button  TEST  &  CALIBRATION  METER(*):  000.0  to  ±199  9  C 

selection  •  l  >  digital  panel  meter  for  monitoring  exposure  calibration  ar 

oiitdiit  unDMti  ~  i  ,  ,  ,  various  switch-selected  (and  probed)  test  points  with 

OUTPUT:  NORMAL  or  COMPLEMENT  data  polarity  selec-  (hp  rprnrrlpr  j 

♦ion.  and  LIN  EAR  or  LOG  code  to  exposure  transfer  charac- ,  ■  •  —  1  ■  .•  •  rr  ,m~ 

teristic  selection  ’  '  "r  r.-T  HIGH-VOLTAGE  SWITCH(*):  enables  high-voltage  powi 

►  *Yi  . _ i _ _ _ j  i _ _ .  •  .  _ 


Dcertiimmi  n.  u  u  ^  t  *  | ■  ...  .  supply  to  be  switched  off  for  safety  during  internal  mai 

RESOLUTION:  Push-button  selection . of; High/, Medium, !  tenance  ■'  ■••’*  -•*>•  . 

or  Low??'  >.  -T,YTv*'  ‘  '  ‘3  , 


evdacmdc.  ,  „  i  __  .  TEST  SELECT  SWITCHES(*):  two  12-position  select 

E’# 10"  ockm9  ype dia  beam  In'-' switches  select  various  test  voltages  and  currents  with 

tensity  adjustment  **•*'■  *"  *  '  the  recorder.  •  Hfr  ■  .u,  *7 


C°NTR°V  °f,f 1 RATE°r  TEST  01096  selec,l0rr  END  FILM  PRIME  MODE  SWITCHC):  enables  CRT  to  be  primed 
and  FILM  ADVANCEmdicators  (disconnected  for  sheet  ,.  mimmize  ima  e  retentjon  .  • 

film  models),  RECORD  push  button  for  initiating  test,  vv  -  v‘ ' 


pattern  or  Image  exposure. 


•  QppY  available  dofti' 
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Appendix  E 


ADDITIVE  NOISE  PRODUCED  BY  FILM  RECORDER 

The  first  nine  steps  of  a  typical  SO-015  duplicate 
(#1102)  were  scanned  with  the  Kitt  Peak  National  Observatory 
(KPNO)  microdensitometer  (10  urn  aperture  at  the  film,  50X 
magnification,  matched  0.25  numerical  aperture  objectives). 
Table  3  lists,  for  each  step,  the  diffuse  density,  D^ , 
measured  with  a  Macbeth  densitometer,  a  number,  NKPN0» 
proportional  to  the  mean  (specular)  transmittance  recorded 
by  the  microdensitometer  and  a  number  N  proportional  to  the 
root  mean  square  fluctuation  in  transmittance.  The  final 

entry,  N.n,  is  the  input  number  to  the  look-up  table  for 
that  step.  From  these  data  we  deduce  that 

log  NKPN0  -  1*12  Dd(SQ_015),  (El) 

and 

loq  NKPN0  =  l-fi8  +  °-85  lc9  Nin'  1  <  Nin-<  110'  (E2) 

Equation  (E2)  somewhat  overestimates  the  measured  values  of 

1°9  f°r  the  first  step  (by  0.05).  This  discrepancy 

apparently  arises  from  the  mismatch  of  the  I UT  and  the 
actual  2415-015  processing,  which  produced  a  somewhat 
greater  slope  than  anticipated  for  small  values  of  N-n>  we 


greater  slope  than  anticipated  for  small  values  of  N.^.  we 
use  an  average  of  the  six  steps  (3  through  8)  to  determine 
a  ,  with  the  aid  of  Eq .  (E2)  and  (E3)  below: 

Nin/cN.  =  (0.4343/od  )  (d (D  t) ) / (d ( log  Nj  ) ) 

in  out 

'NKPNo/on  >  (d(log  NKPNO)/d(1°^  Nin))#  (E3) 

K  PN  U 

Eq.  (E3)  makes  use  of  the  fact  that  AD  /o_.  is 

^  out'  D 

out 

independent  of  the  instrument  used  to  measure  density  and 
density  fluctuations  {whether  diffuse  or  specular),  and  that 

NKPNO  *s  proportional  to  the  (specular)  transmission  of  the 
wedge.  Thus,  with  Eq .  (E2)  and  (E3)  and  the  entries  of 

Table  13,  we  generate  the  following  values  of  a ^  for  steps  3 
through  8: 


TABLF  n 


Dicomed  Noise  (e^) 

versus 

iuuut  Number  (iv 

Step  # 

D 

N. 
l  n 

nkpno 

cNkpno 

a  N 

3 

0.65 

34 

12.  59 

3.18 

4 

0.50 

51 

17.88 

3.  36 

5 

0.41 

68 

23.64 

3.  39 

6 

0.  33 

85 

30.63 

3.  27 

7 

0.  26 

102 

35.89 

3.  34 

8 

0.  21 

119 

40.27 

3.47 

From  these 

entries  we 

find 

<v  - 3- 

33,  <(oN  - 

<0  >) 2S  ^  0.10 

N 

P'-M) 

We  conclude 

that  an 

additive,  uniform, 

uncor  rel  -i  i.od 

:ootr ibution 

to  the 

input  noise  is 

an  excellent 

’!vtn  •i>*rization  of  the  Dicomed  noise.  Its  rms  value  is 
3.33  on  an  8-bit  scale. 

The  actual  input  SNR  for  the  hardcopy  is  accordingly 
reduced  as  indicated  in  Table  14. 
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TABLE  14 

Input  SNR  £or  Hardcopy 


Noise  Level  SNRin  ~  Normal  SNRin  ~  Hardcopy 

0  200  75.1 

1  100  60.0 

2  50  41.6 

3  25  23.7 


4 


12.5 


12.3 


Appendix  F 

PHOTOGRAPHIC  PROCESSING  FOR  HARDCOPY 

The  maximum  densities  measured  from  the  Precision  Photo 
originals  are  listed  in  Table  15. 


TABLE  15 

Maximum  Densities  in  Original  Photo,  by  Scene 


Scene 

Density 

Scene 

Density 

01 

2.1 

81 

1.8 

11 

1.8 

91 

h-* 

• 

00 

41 

1.7 

12 

1.9 

51 

2.1 

42 

1.7 

61 

1.7 

52 

2.0 

We  intended  that  the  maximum  density  of  our  final 

hardcopy  (DQUt  )  be  approximately  1.8,  this  being  the 
max 

average  value  of  the  Precision  Photo  originals  which  we 
wished  to  emulate.  Experience  in  2415/SO-015  contact 
printing  and  Dektol  processing  warned  us  that  the  post- 
Dicomed  hardcopy  originating  from  the  darkroom  at  USC  would 
have  to  be  no  higher  in  density  than  around  1.6.  We  also 
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learned  that  we  could  successfully  "contrast  stretch" 
maximum  densities  between  1.3  and  1.6  to  our  goal  of  1.8  by 
careful  Dektol  processing;  however,  Dicomed  playbacks  darker 
than  1.6  inevitably  resulted  in  higher  SO-015  densities. 
Calibration  was  achieved  by  contact-printing  (30  second 
exposure)  a  Kodak  stepwedge  onto  SO-015,  which  was  then 
developed  in  Dektol  at  68°F  for  times  ranging  from  2  to  7 
minutes.  Calibration  data  are  illustrated  in  Figure  26.  We 
next  drew  curves  of  input  density  against  development  time 
to  yield  constant  output  density,  and  calculated  a 
correction  for  the  contact  exposure  time  in  order  to  shift 
the  actual  maximum  2415  density  to  a  suitable  value  in  the 
abscissa.  This  latter  constraint  was  designed  to  give 
development  times  and  shifts  of  relative  input  density  which 
would  produce,  as  near  as  possible,  linear  2415/SO-015  tonal 
transfer  curves. 

Table  16  lists  the  results  of  the  final  maximum  hardcopy 
stepwedge  densities.  Our  early  attempts  to  accommodate  the 
variability  of  the  Dicomed-2415  output  are  reflected  in  the 
relatively  high  values  of  mean  and  standard  deviation  within 
each  scene  for  blur  levels  2  and  0.  (Although  stepwedge 
data  are  not  available  for  blur  level  2  images,  we  were  able 
to  estimate  the  maximum  density  from  comparisons  with  trends 
in  other  blur  levels.)  Contributing  to  the  improved 
processing  was  our  realization  of  the  need  to  adjust  the 
exposure  time  during  contact  printing  for  each  transparency. 


89 


L_ _ 1 _ 1 _ l _ 1 _ 1 

0 

0.5 

1.0  1.5  2.0  2. 

°2l.l5  ~ 

?«*0 

160 

80  0 

J 

"*  NOUt 

I  1  1 

HIN 

♦  0.6 

♦  1.2  ♦l.fi  ♦2.«t 

2. *4 

t.f 

0.8  0.0 

LOG 


Figure  26:  Characteristic  curves  for  S0-015  film  developed 
in  Dektol. 


The  derivative  of  the  characteristic  of  each  curve  is  also 
illustrated.  The  approximate  (ideal)  relationship  between 
D  and  the  quantities  of  Figure  17  is  illustrated  by  the 
additional  three  scales. 


instead  of  working  with  batches  at  constant  exposures, 
thereby  allowing  the  2415  variability  to  be  compensated  by 
the  printing  and  development.  The  densities  of  Table  16  are 
average  values  for  the  set  of  five  duplicates  that  were 
produced;  the  variability  over  these  copies  was  similar  to 
that  shown  within  the  five  noise  levels. 
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TABLE  16 


912 

1.50 

1.70 

913 

1.74 

1.67 

914 

1.85 

1.82 

120 

2.03 

1.89 

121 

1.60 

1.66 

122 

2.96 

3.00 

123 

1.75 

1.77 

124 

1.88 

1.88 

420 

1.82 

1.90 

421 

1.86 

1.89 

422 

1.65 

2.25 

423 

1.92 

1.91 

424 

1.82 

1.82 

520 

1.89 

1.93 

521 

1.76 

1.74 

522 

1.10 

1.10 

523 

1.78 

1.80 

524 

1.78 

1.76 

1.80 

1.85 

1.80 

1.70 

1.76 

1.69 

1.84 

1.81 

1.81 

1.74 

1.78 

1.75 

1.65 

1.64 

1.71 

2.47 

2.72 

2.28 

1.73 

1.75 

1.70 

1.88 

1.86 

1.88 

1.96 

1.97 

1.78 

1.91 

1.94 

1.89 

2.55 

2.80 

2.35 

2.01 

1.96 

1.93 

1.87 

1.79 

1.82 

2.32 

1.97 

1.60 

1.79 

1.73 

1.76 

1.20 

1.20 

0.90 

1.73 

1.77 

1.75 

1.75 

1.74 

1.75 

1.73 

0.14 

1.8 

1.71 

0.04 

1.83 

0.02 

1.84 

0.12 

1.65 

0.04 

2.69 

0.32 

1.8 

1.74 

0.03 

1.88 

0.01 

1.89 

0.08 

1.90 

0.03 

2.32 

0.43 

1.95 

1.95 

0.04 

1.80 

0.02 

1.94 

0.26 

1.76 

0.02 

1.10 

0.12 

1.8 

1.77 

0.03 

1.76 

U.01 

We  decided  at  the  time  of  processing  to  stretch  the 
stepwedge  dynamic  range  of  scenes  41,  42,  and  61  beyond  1.8. 
This  was  prompted  by  their  broader  histograms  (particularly 
scene  41)  as  evidenced  in  the  lighter  appearance  of  the 
scenes. 

Over  the  range  of  input  densities  in  2415  between  0.0  and 
1.5,  the  three  calibration  curves  of  Figure  1  are  accurately 
described  as 

°015  a  a  +  b  D2415  +  c  D22415  (pl) 

where  a,  b,  and  c  are  taken  from  Table  17  according  to  the 
appropriate  development  time,  tD. 

The  three  additional  scales  parallel  to  the  D2415  scale 
in  Figure  26  illustrate  approximately  the  relations  among 
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TABLE  17 


Coefficients  for  Each  Development  Time 


tD(min) 

a 

b 

c 

3 

0.10 

0.16 

0.51 

4 

0.  11 

0.  20 

0 . 66 

5 

0. 13 

0.25 

0.74 

the  values  of  D.  the  density  of  the  Precision  Photo 
original s,  log  Nin,  the  log  of  the  corresponding  digital 
value  recorded  by  the  microdensitometer  (truncated  to  8 
bits) ,  and  Nout*  the  number  presented  to  the  Dicomed 
recorder.  The  weakest  approximation  is  that  between  N 

out 

and  D24i5*  The  N0ut  ran9e  was  randomly  stretched  and 
compressed,  with  respect  to  the  D24i5  scale  in  the  hand¬ 
processing  of  the  2415  sheets  at  USC.  The  reason  for  this 
variability  is  apparent  from  a  study  of  Figure  27,  taken 
from  a  Kodak  data  release.  The  characteristic  curves  of 
2415  film  (then  called  SO-115)  are  extr2mely  sensitive  to 
small  changes  in  development  conditions.  These  curves 
illustrate  changes  with  development  time.  Temperature 
changes  would  be  equally  dramatic,  as  would  agitation 
conditions. 

We  asked  USC  to  work  toward  a  maximum  2415  density  of  1.6 
in  the  step  of  the  wedge  corresponding  to  N  *  0,  but  to 
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LOG  EXPOSURE 
(METER-CANDLE-SECONDS) 

Figure  27:  Characteristic  curves  for  Kodak  Technical  Pan 
2415  Developed  in  Kodak  D-19  at  68°F  developed 
in  Kodak  D-19  at  68°  F 


err  when  necessary  (it  was)  toward  underdevelopment.  As 
Figure  26  illustrates,  such  errors  could  be  compensated  by 
reducing  the  exposure  of  the  015  copies  and  increasing  their 
development  time.  Both  actions  would  increase  the  slope  of 

the  D015  versus  D2415  curve  for  the  actual  2415  record, 
though  generally  not  enough  to  compensate  for  the  reduced 
contrast  in  the  highlights  attendant  to  underdevelopment  of 
the  2415  sheets. 


